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FOREWORD 


This report covers work on one phase of the photoheliograph 
development task, NASA Code 945-84-00-01-00, for the period 
November 1967 through June 1968. The photoheliograph has 
been proposed to NASA for the Apollo telescope mount (ATM) by- 
Caltech, with Professor Harold Zirin as the principal investi- 
gator and Dr. Robert Howard of Mt„ Wilson and Palomar 
Observatories the co -investigator (see TM 33-369, November 
1967). The objective of the investigation is to obtain high reso- 
lution cinematographs in white light near ultraviolet and narrow 
band hydrogen alpha. Because of the ATM program uncertain- 
ties, emphasis has been placed on areas of technology that are 
somewhat mission-independent, but the ATM spacecraft has 
been used to establish design constraints. 
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PHOTOHELIOGRAPH 
PRIMARY MIRROR DEVELOPMENT 

In a near -diffraction limited system, the ultimate success of the system 
depends on the proper choice and fabrication of materials for the optical' ele- 
ments. The primary mirror functional requirements are listed in Table 1. 
These requirements necessitate use of a primary mirror material of extremely 
high optical quality and mechanical dimensional stability and low thermal dis- 
tortion under solar heating conditions during service in orbit. 

In order to select an optimum material to meet these requirements, an 
index, called the Thermal Distortion Index (TDI) was used for evaluating the 
thermal-physical merits of various candidate optical primary mirror mate- 
rials. This index is independent of geometric configuration and is a function 
only of the physical properties of the materials. The TDI is a qualitative 
measure of the ability of a material to resist a change in its dimensions with a 
temperature gradient applied across the material. The TDI is the ratio of the 
coefficient of thermal expansion to the thermal diffusivity or 


TDI 



where 


a - coefficient of Linear thermal expansion 
K = thermal conductivity 
"6^= specific heat 
p = density 

A low numerical number is desired for the TDI because materials having 
a low TDI have a lesser tendency to support temperature gradients and will 
stabilize quicker with Less optical distortion, Also a high elastic modulus to 
density ratio (E/p) makes possible a stiffer structure for a given weight. A 
low TDI ratio coupled with a low elastic Limit and a low damping capacity would 
point toward problems during launch. Pertinent physical properties and TDI's 
of a number of candidate materials are given in order of decreasing index in 
Table 2. 
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Table 1. Photoheliograph, Primary Mirror Functional Requirements 

The Photoheliograph primary mirror design and fabrication shall 
assure the following: 

1. Systematic wavefront accuracy of the total optical system 
is l/8th wave or better over a two - - arc -minute half-field, 
and random, small-scale errors are l/50thwave or less 
during operation of the telescope. 

2. Thermal distortion does not degrade the systematic error 
by more than l/20th wave during operation of the 
telescope. 

3. Vibration, shock, acoustic, and thermal cycling environ- 
ments do not permanently degrade the systematic error to 
more than that specified in Requirement 1. 

4. During telescope operation, the primary mirror can be 
rigidly supported within ±0. 004 in. of a predetermined 
"aligned” position with the optical axis tilted no more than 
0.5 arc-minutes from the predetermined "aligned" 
condition. 

5. Operation in a l-g_ gravity field will distort the mirror 
surface no more than 1/40 wave from the surface figure 
expected in a 0-j* field. 
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Table 2. Comparison of Thermal Performance for Various Candidate Primary Mirror Materials 


Material 

a 

(X0 -6 /°F) 

P 

(Lb /in?) 

C 

P 

(Btu/lb- °F) 

K 

(Btu/Hr-Ft- °F) 

TDI 

(10" 10 Hr-Ft/°F-in?) 

E 

(10 6 lb /in. 2 .) 

E/p 

(10 6 in. ) 

300 Series Stainless Steel 

9.2 

0.29 

0. 12 

9.4 

340 

28 

97 

Titanium (Com'l. pure) 

5. 0 

0. 16 

0. 13 

9.8 

106 

16.5 

103 

Nickel (Com'l. pure} 

7. 4 

0. 32 

0. 11 

45. 0 

58 

30 

94 

Fused Silica ( 7940 )< 1 ) 

0. 28 

0. 079 

0. 17 

0. 77 

49 

10.6 

134 

Brass (70 Cu, 30 Zn) 

11. 1 

0.31 

0. 09 

70. 0 

44 

16 

52 

Aluminum (5052) 

13. 2 

0. 097 

0. 23 

80. 0 

37 

10.2 

105 

Magnesium (ZK60) 

14.5 

0. 066 

0. 24 

69. 0 

33 

6. 5 

98 

Tantalum (Com'l. pure) 

3. 6 

0.60 

0. 036 

31. 0 

25 

27 

45 

Invar 

0.5 

0.29 

0. 12 

7.8 

22 

21 

72 

Beryllium 

6.4 

0. 066 

0. 45 

104. 0 

18 

42 

635 

CER-VIT C-101^ 

0. 08 

0. 090 

0.22 

0. 97 

16 

13.4 

149 

Copper (Com'l. pure) 

9.3 

0. 32 

0. 092 

225. 0 

12 

16 

50 

Silver (Com'l. pure) 

10.6 

0. 38 

0. 056 

242. 0 

9.3 

10.5 

28 

U.L.E. Fused Silica (7971)t 2 > 

0. 03 

0. 079 

0. 18 

0. 75 

5. 7 

9.8 

124 

Tungsten (Com'l. pure) 

2. 4 

0. 70 

0. 034 

116. 0 

4.9 

50 

71 

Super Invar 

0 . 1 

0. 29 

0. 12 

8. 0 

4. 4 

21 

72 

Graphite (Molded) 

1.4 

0. 064 

0. 20 

70. 0 

2.5 

0. 7 

11 

^ Corning Glass Works amorphous fused silica glass. 

(^Corning Glass Works amorphous ultra Low expansion (U.L.E.) fused silica glass. 

^Owens Illinois crystallized glass. 

NOTE: Table compiled by R. A. Happe 
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A number of materials can be eliminated on the basis of the following 
rationale. Stainless steel {300 series) and titanium have very high TDI. 

Nickel and brass have moderately high TDI and are not competitive with copper. 
Magnesium does have a low TDI, but is not competitive with aluminum with 
regard to figuring coatings, corrosion resistance and fabricability such as 
brazing. Tantalum is expensive and has a low E/p ratio. Copper has a low 
E/p ratio, is too soft to figure, and compatible hard figuring coatings are not 
available. Silver has a very low E/p and appears to have no particular advan- 
tages over copper. Tungsten has a low TDI but is expensive, very difficult to 
machine, and material availability for the size of the primary mirror is ques- 
tionable. Graphite has a very low E/p ratio and would be almost impossible to 
figure to the near diffraction limited tolerances required. 

The surviving candidates are listed below. 

Fused Silica (7940) 

Aluminum (5 052) 

Invar 

Beryllium 

CER-VIT C-101 

U. Li. E. Fused Silica (7971) 

Super Invar 

A static (equilibrium) thermal distortion computation was made for each 
of the candidate materials in a flat mirror blank configuration. The configura- 
tion was a one-in. -thick, 31. 5-in. -dia. blank, simply supported around the 
entire periphery. The back face was held at 70 °F; the front face with a solar 

absorptance ( a ) of 0, 14 was sun illuminated. The maximum equilibrium dis- 
s 

placement of the front surface was computed for each material. The materials 
with the minimum steady-state deflections were Super Invar, U.L. E, fused 
silica (7971), and CER-VIT C-101. 

Concurrently with the above analysis, a study was made to determine the 
effect of primary mirror configuration on the maximum steady- state deflection. 
Figure 1 compares the displacement of a fused silica flat mirror with configu- 
rations thicker at the outside, thinner at the center, and vice versa. The con- 
figuration corresponding to the bottom curve of Figure 1 was chosen for further 
analysis. 
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Figure 1. Effect of Primary Mirror Design on Thermal Distortion, F/49 Gregorian Telescope 
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Further investigation revealed that although Super Invar was developed in 
the early 1930's, little was known about its metallurgical and fabrication char- 
acteristics. It was felt that a development program would have to be accom- 
plished before Super Invar could be considered as a candidate for the primary 
mirror. This development program would have to include as a minimum, the 
following: 

1. Investigate the effects of impurities on the coefficient of linear 
thermal expansion (a) and the inclusion count of annealed, wrought 
material. 

2. Develop optimum heat treatment for minimum a and optimum 
dimensional stability. 

3. Develop figuring techniques to the stringent optical tolerances 
required for the primary mirror. 

The cost of this program was estimated to be in the neighborhood of $225K 
which could not be justified within the constraints of the program. 

The two materials left for consideration as candidates for the primary 
mirror were U. L.E. fused silica and CER-VIT G-101, Throughout the pre- 
liminary development phase of the primary mirror, serious consideration had 
been given to an active internal cooling system to minimize the mirror optical 
distortions and to reduce the heat flux from the primary mirror system such 
that the overall telescope heat transfer to the ATM spar was less than 
0.6 Btu/hr. Since 50 watts are calculated to be absorbed by the primary 
mirror, the need arose for active internal cooling to maintain its equilibrium 
temperature at approximately the same temperature as the ATM spar. Instal- 
ling coils on the back surface of the mirror necessitated use of a very low 
expansion material such as Super Invar or possibly U.Li„E. fused silica itself. 
The problem of the attachment of the coils had no clear solution and the result- 
ant thermal and mechanical distortions were much higher than required for 
internal cooling the mirror on its neutral axis. Radiative cooling required that 
the back surface of the mirror viewed a refrigerated surface (heat sink) with an 
emissivity of 0. 9 held at -40 °F. The heat sink surface would require advanced 
development work to attain a refrigeration cycle capable of operating reliably 
for a long time in a zero g field. 

A design evolved for locating the cooling coils on the neutral axis of a • 
65- cm primary mirror corresponding to configuration No. 1 of Figure 2. The 
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CONFIGURATION No. 3 


Figure 2. Primary Mirror Configuration Identifications 
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fabrication concept involved using amorphous U.L.E. fused silica, installing 
the cooling channel grooves in the top piece, and then fusing the top and bottom 
plates together, 

CER-VIT C-101 is a crystalline glass and cannot be fused to itself and 
still maintain its desired physical properties after fusing. Lengthy discussions 
were held with the Owens-Illinois technical personnel for practical solutions to 
an active internally cooled CER-VIT C-101 mirror. A telecon with 
Dr. Simmons on 21 March, 1968 revealed that Owens-Illinois chose not to bid 
on manufacturing a primary mirror blank from CER-VIT C-101 with internal 
cooling coils. Because of limited time and manpower other techniques for 
active cooling a CER-VIT C-101 primary mirror blank other than internally 
could not be pursued. 

The effects of gravity sag on configuration No. 1 for U.L.E, fused silica 
were checked and it was found that the 1/2 optical path difference (1/2 O.P. D. ) 
was a nineteenth of a wave (see Table 3). The structural analysis is described 
in Appendix A. The O.P. D. is actually twice this or \/9. 5, which is clearly 
unacceptable (see Item 5, Table 1, which requires a minimum O.P.D. of X/40 
due to gravity g sag). The 1/2 O. P. D. is illustrated schematically in Fig- 
ure 3 and is self-explanatory. The 1/2 O.P.D. calculation is described in 
Appendix B. A sample 1/2 O.P.D. calculation is given in Table 4. 

The effect of g sag on four other configurations (see Figure 2) was calcu- 
lated for U.L.E. fused silica and the results tabulated in Table 3. Two satis- 
factory shapes were configurations 2 and 5. Although configuration 5 gave a 
slightly better 1/2 O.P.D. , its parabolic back face would be more difficult to 
fabricate than the double tapered back of configuration 2. Therefore configura- 
tion 2 was considered for further analysis with U.L.E. fused silica. 

Recent work on beryllium dimensional stability from Battelle Memorial 
Institute indicated that grades with high BeO contents (4-5%) vs 2% for the 
former S200D grade and a smaller grain size, 10 microns average with no 
single grain exceeding 15 microns, increased the precision elastic limit, 

P. E. L, (the stress required. to induce a permanent strain of lpin. /in. ) from 
8000 psi to as high as 40, 000 psi. This increased dimensional stability led to 
speculation of its candidacy again as a primary mirror blank material. Gravity 
sag calculations for configurations 2 and 5 revealed that the effect of j* sag was 
insignificant (see Table 3), 
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Material 


Primary Mirror 
Configuration 


1/2 Optical Path 
Difference 
{1/2 OPD) 


Focus Change 
(Inches) 


U.L. E. Fused Silica 
U. L. E. Fused Silica 
U. L. E. Fused Silica 
U.L. E. Fused Silica 
U. L. E. Fused Silica 
Beryllium 


1 

2 

3 

4 

5 
5 


X 

19 


X 

185 

X 

48 


X 

44 


X 

25 0 


X 

1000 


0. 0051 
0. 0030 
0. 0031 
0. 0023 
0. 0021 
0. 0005 


Beryllium 


2 


X 

909 


0. 0006 











Table 4. Calculation of 1/2 O.P.D. for ULE Fused Silica, Configuration No . 1, 1G Load in -Y Direction 


X 

X 2 

M 

>4 

11 

Ay 

(10“ 6 in. ) 

Ay - Ay 0 
Ay =0 @ X=0 
(10 -6 in. ) 

y" 

= y + (Ay -Ay Q ) 

y' 

- (y' - y"> 10 ^ in. ) 

0 

0 

0 

+0.088 

0 

0 

0 

0 

2 

4 

0.010160000 

-0.286 

-0.374 

0.010159626 

0.010159469 

-0. 157 

4.32 

18.66240 

0.047402496 

-2.509 

-2.597 

0.047399899 

0.047400018 

+0.119 

6.04 

36.48160 

0.092663265 

-5.238 

-5.326 

0.092657939 

0.092658421 

+0.482 

8.07 

65.12490 

0. 165417247 

-9.381 

-9.469 

0.165407778 

0. 165408600 

+0.822 

9.79 

95.84410 

0.243444016 

-13.532 

-13.620 

0.243430396 

0.243431289 

+0.893 

11.82 

139.71240 

0.354869499 

-18.953 

- 19-041 

0.354850458 

0.354850947 

+0.489 

12.79 

163.58410 

0.415503618 

-21.634 

-21.722 

0.415481896 

0.415481896 

0 


o 

t 


/ _ _ 163.58410 _ 40.896025 _ qR 4 , n , 4171 

1 (@X=12.79 in.) " 0 . 415481896 " 70 


/ _ X 2 _ 

7 ~ 4f' ~ 4 (98.43034171) 


= 0.002539867236X 


Af = f' - f = 0.0051 in. 

ATX = 5000 A, X = 20 x 10 ^ in. (20 micro inches), or 1 microinch (1 |i in.) = X/20 
Maximum deviation above (y -y " ) is 1.05 [j. in. 


1 p. m. 
_X_ 
20 


. 1 • 05 F ... x = 19 . 05 , the 1/2 O.P.D. is .*. 

IT 
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Figure 3. Definition of 1/2 OPD 
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A detailed thermal analysis was made of configuration 1 in U.L.E. fused 
silica before the computer results of the g-sag distortions of the other configu- 
rations were made available. These results are tabulated in Table 5. The 
thermal analysis is described in Appendix C, The coolant used was 80% 
methanol, 20% water which entered the mirror at 70 °F at 16 psia. It can be 
seen that turbulent flow gives slightly better results than laminar flow and that 
the coolant fluid flow direction is better from the mirror inside to the mirror 
outside. Because of the unknown possible problems of vibration induced by 
turbulent flow it was decided to choose the conservative {lower flow rate) 
laminar flow, with the flow from inside to outside. It was recognized that tur- 
bulent flow gave a TF radial temperature rise from the mirror inside to out- 
side where laminar flow gave a 3°F temperature rise. The 1/2 O. P. D. for 
U.L.E. fused silica corresponding to the thermal condition chosen above for 
configuration 2 was a satisfactory \/55. The local deflection of the U.L.E. 
fused silica mirror reflecting surface due to cooling fluid pressure (16 psi) 
corresponded to a 1/2 O.P.D. of X/500, and the deflection of the reflecting 
surface from the 70 °F laminar flow temperature distribution (thermal ripple) 
corresponded to a 1/2 O.P.D, X/2000, both insignificant. 

Because of the known fabrication difficulties of installing cooling channels 
internally in a beryllium mirror blank, and because of the limited time avail- 
able, it was decided to compare a back face radiatively cooled beryllium mir- 
ror blank with a U. L. E. fused silica mirror blank both of configuration 2, the 
back face double taper. The results (Table 5) show that for a temperature gra- 
dient on the mirror back, 3°F rise from inside to outside, the 1/2 O. P. D. for 
U.L.E. fused silica was X/20 and the 1/2 O.P.D. for beryllium \/ 4. Moving 
the coolant to the back of the mirror not only loses the advantages of close 
proximity to the heated surface, but does not avoid problems. The beryllium 
could not be considered for further analysis because of the poor optical per- 
formance (even X/ 20 is marginal). Another problem with beryllium is that a 
nickel-phosphorous Kanigen coating would be required for figuring. This coat- 
ing presents a poorer optical surface than glass and is insufficient for use in 
the near ultraviolet as UV- scattering results. 

Table 6 lists the’re suits of the transient thermal distortion analysis for 
U.L.E. fused silica, configuration 2, with laminar flow from mirror inside to 
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Table 5. Distortic 


Material 

Primary 

Mirror 

Configuration 

Thermal Condition 

l/2 Optical Path 
Difference 
(1/2 OPD) 

Focus Change 
(Inches) 

U. L. E. Fused 
Silica 

1 

Laminar flow, uniform 
70 °F coolant temperature 

X 

57 

0.0008 

U. L. E* Fused 
Silica 

1 

Turbulent flow, uniform 
70 °F coolant temperature 

X 

78 

0.0006 

U. L. E. Fused 
Silica 

1 

Turbulent flow from 
outside to inside 

X 

72 

0.0006 

U.L.E. Fused 
Silica 

1 

Turbulent flow from 
inside to outside 

X 

83 

0.0007 

U. L. E. Fused 
Silica 

1 

Laminar flow from 
outside to inside 

X 

49 

0.0007 

U.L.E. Fused 
Silica 

1 

Laminar flow from 
inside to outside 

X 

Vf 

0.0010 

U. L. E. Fused 
Silica 

2 

Temperature gradient on 
mirror back, 3°F rise 
from inside to outside 

X 

20 

0.0025 

Beryllium 

2 

Temperature gradient on 
mirror back, 3°F rise 
from inside to outside 

X 

4 

0.0417 

U. L. E. Fused 
Silica 

2 

Laminar flow from 
inside to outside 

X 

55 

0.0007 
















Table 6. Distortion Analysis - Transient Thermal 


Material 

Primary 

Mirror 

Configuration 

Thermal 

Condition 

Minutes of 
Direct Solar 
Illumination 

1/2 Optical Path 
Difference 
(1/2 OPD) 

Focus Change 
(Inche s) 

U.L.E. Fused 
Silica 

2 

Laminar flow 
from inside 
to outside 

12 

X 

49 

0. 0011 

U.Li.E. Fused 
Silica 

2 

Laminar flow 
from inside 
to outside 

24 

X 

59 

0. 0010 

U.LuE. Fused 
Silica 

2 

Laminar flow ' 
from inside 
to outside 

36 

X 

61 

0. 0009 
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outside. It can be seen that the optical performance is increasing from 12 to 
36 minutes solar illumination -time from passing from earth shade to direct 
solar illumination; the optical performance for each time-period is acceptable. 
The mirror starts the sun portion of the orbit at 70 °F. The coolant fluid starts 
at this temperature by means of a fluid temperature controlled loop as shown in 
the schematic in Figure 4. The schematic is correct except for the fact that 
the fluid should enter at the innermost channel at 70 °F, 

The U.L.E. fused silica mirror blank containing the integral cooling 
channels is fabricated in the following way. An approximate spiral channel 
would be sandblasted into the forward face of the bottom plate of the mirror, 
the channels fire polished, and the top (approximately one-half inch thick) and 
bottom plates fused together. The procedure would be to heat slowly to 
1000 °C, then quickly to 1600 °C for fusing, and then sag over a refractory mold 
to general shape at just under 1600°C. The sagging process might introduce 
sufficient change of the front face contour from that desired such that subse- 
quent optical grinding of the front face could decrease the 0.5-in. nominal 
thickness of the mirror front plate by as much as 0. 15 in. This will result in 
a negligible change of 1/2 O.P.D. (X./56 vs \/55 1/2 O.P.D. for a uniform 
0. 50 inch front plate thickness). The mirror blank would then be rough ground, 
and annealed at 1000°C. The U.L.E. fused silica mirror blank fabricator, 
Corning Glass Works, will guarantee 90% fusion with 100% fusion around the 
outside periphery. It is planned to let the purchase obligation for one primary 
mirror blank flight verification unit with integral cooling channels. 

One disadvantage of U.L.E. fused silica is that high internal stress 
causes some concern optically and mechanically. The internal stress arises 
from inability to completely anneal U. L, E. fused silica. These stresses cause 
optical birefringence. The magnitude of the measured birefringence, there- 
fore, is a measure of the internal stresses. Corning Glass Works will only 
guarantee a value of 50 millimicrons per centimeter of transmitted light path. 
One other disadvantage is that it is a fairly new material and the fabrication, 
testing, and figuring techniques probably have not been optimized yet. 

CER-VIT C-101, the Owens-Illinois non-fusable crystallized glass has a 
minimum guaranteed value about 1/5 the birefringence that the Corning product 
has. It has been on the market about two years longer than U.L.E, fused 
silica and the optical figuring people have gained much more experience in 
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Figure 4. Fluid Temperature Controlled Primary Mirror Schematic 






750-7 


working with this material. With a premium select grade of CER-VIT C-101, 

a coefficient of thermal expansion can be obtained virtually identical to that of 

-6 

U.L. E. fused silica, i. e. , 0.03 x 10 /°F. Because of these desirable quali- 

ties, engineering discussions shall be reopened with Owens-Illinois concerning 
the problem of thermal coupling to the premium select CER-VIT C-101 mirror 
blank, if it appears Corning is unable to fabricate a U.L.E. fused silica pri- 
mary mirror blank. The design shall be specifically developed for a mirror 
blank without internal cooling. 
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APPENDIX A 

STRUCTURAL ANALYSIS OF THE PHOTOHELIOGRAPH 
PRIMARY MIRROR 

INTRODUCTION 

The objective of the report is to summarize the analytical effort to date 
on the structural analyses of the primary mirror. The analyses consist of the 
thermal and gravitational distortion of various mirror configurations and 
materials . 

First, a brief discussion of the theory and assumptions are presented; 
second, the effort to verify the analysis method by solution of example prob- 
lems is presented; and third, the results of all analysis to date are summarized 
in Table 7 . 

THEORETICAL DEVELOPMENT 

The analyses were performed by use of a digital computer program titled 
n ELAS"{l). The theoretical description to follow includes only that pertinent to 
the mirror. 

A finite element approach is used to characterize the structure of the 
mirror. The mirror is discretized into small finite structural elements; within 
each element the material properties, loading conditions and displacement 
characterization are assumed constant. A good approximation to the true 
answer is obtained by an increase of the number of elements in the mirror 
where the displacement gradients are high. 

The geometric and load axisymmetry of the mirror reduced the analysis 
to a two-dimensional problem. A triangular torus and quadrilateral torus 
finite elements are used to represent the mirror. 


^Utku, S., and Akyz, F.A. , ”ELAS - A General-Purpose Computer Program 
for the Equilibrium Problems of Linear Structures”, JPL TR 32-1240, VoLl, 
Feb. 1, 1968. 
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Table 7. Summary of the Computer Analyses on the Photoheliograph Primary Mirror 


Hun 

No, 

Conf . 

Material 

Conditions 

Maximum 
Effective 
Stress 
(P.S.I. ) 

Maximum 

Stress 

Node 

Location 

Maximum 
Defl in 
V Direction. 
Inches x 10“® 

Maximum 
Defl. in 
XDirection 
Inches x 10“ 6 

1 

1 

ULE Fused Silica 

Laminar flow, uniform 70°F coolant temperature 

0. 55 

54 

-3.489 

1.867 

2 

1 

TILE Fused Silica 

1 G load, mirror support down 

4. 62 

30 

-22.372 

2.478 

3 

1 

ULE Fused Silica 

1 G load, mirror support up 

3. 84 

21 

22.372 

-2.478 

4 

1 

ULE Fused Silica 

Turbulent flow, uniform 70° coolant temperature 

0. 41 

100 

-2.768 

0.929 

5 

1 

ULE Fused Silica 

Turbulent flow from outside to inside 

0. 52 

65 

-2.422 

1.094 

6 

1 

ULE Fused Silica 

Laminar flow from outside to inside 

0.78 

52 

-2.499 

2 366 

7 

1 

ULE Fused Silica 

Laminar flow from inside to outside 

0. 50 

86 

-4.294 

2 454 

8 

1 

ULE Fused Silica 

Turbulent flow from inside to outside 

0.35 

86 

-3.051 

1.126 

9 

2 

ULE Fused Silica 

1 G load 

18. 05 

47 

-12.938 

1 782 

10 

3 

ULE Fused Silica 

1 G load 

10.45 

47 

-13.485 

1.897 

11 

4 

ULE Fused Silica 

1 G load 

29. 05 

47 

-10.220 

1.723 

12 

5 

ULE Fused Silica 

I G load 

21.52 

47 

-9.640 

1 534 

13 

5 

Beryllium 

1 G load 

17.78 

47 

-2.046 

0 323 

14 

2 

B eryllium 

1 G load 

14.91 

47 

-2.715 

0.375 

15 

2 

ULE Fused Silica 

1 G load, the mirror was supported at the outside 
edge 

8. 44 

261 

-4 233 

-0.671 

16 

2 

ULE Fused Silica 

1 G load, the mirror cooling channels were removed 

18. 51 

47 

-13.123 

1.786 

17 

2 

Beryllium 

1 G load, the mirror cooling channels were removed 

15. 30 

47 

-2,757 

0.376 

18 

2 

ULE Fused Silica 

Temperature gradient on mirror back, 3°F rise from 
inside to outside 

1. 50 

72 

-10,827 

4.293 

19 

2 


Void, incorrect temperature instruction chosen 





20 

2 

B eryllium 

Same conditions as Run No. 18 

393. 50 

241 

- 190.969 

143.129 

21 

2 

ULE Fused Silica 

Laminar flow from inside to outside 

0. 53 

212 

-3.011 

2.341 

22 

2 

ULE Fused Silica 

Laminar flow from inside to outside 12 min. sun 
illumination 

0. 55 

54 

-4.824 

1.938 

23 

2 

ULE Fused Silica 

Laminar flow from inside to outside 24 min. sun 
illumination 

0. 53 

191 

-4.493 

2.213 

24 

2 

ULE Fused Silica 

Laminar flow from inside to outside 36 min. sun 
illumination 

0. 56 

212 

-4.113 

2.289 

25 

2 

CER-V1T 

Same conditions as Run No. 18 

1.92 

71 

-9.777 

3.850 

26 

2 

ULE Fused Silica 

Mirror 15 mm. thinner at center line duo to 
fabrication 

0. 55 

212 

-2.784 

2.292 

27* 








zs» 








29* 








30* 








31 

2 

ULE Fused Silica 

Local deflection of mirror reflecting surface due 
to cooling channel presfluro (Id psi) 

32.85 

77 

0.282 

-0. 171 

32 

2 

ULE Fused Silica 

Local defloction of mirror reflecting surface from 
temperature distribution of 70°F laminar flow 
(thermal ripple) 

1 93 

34 

-0.0955 

0. 136 

*Still in p 

rogress 


750- 















Cross Section of a Mirror 


Half of a Triangular Torus 


Each node is represented by two translational degrees of freedom and the 
relationship between nodal displacements and forces for each element is 
derived. The relationships for all elements of the mirror are assembled, the 
boundary conditions are imposed, the loads are applied, and the equations are 
solved for the deflections and stresses. 

The stiffness finite element approach is used where a displacement func- 
tion for each element is defined to match the nodal displacements and to insure 
displacement continuity between elements. From the displacement function 

and material characteristics the following relation between the nodal displace- 

til 

ments and forces are obtained for each finite element. For the i — element. 



where 


. = sti 

Mi- 

If I,- 


stiffness matrix for i— element 


nodal displacement for i — element 


nodal forces for i ~ element 


The equations of equilibrium for the entire mirror are obtained by generation 
of the total stiffness matrix by the addition of all elemental stiffness matrices, 
elimination of rows and columns of stiffness matrices corresponding to 
restraints, and the application of either thermal or gravitation forces. The 
equation is fit 

MM -I* 
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75 0-7 


where 



= stiffness matrix of total structure 
= unknown nodal displacements 



applied force 


The equation is solved for the nodal displacements by 


W H' 1 i F l 

The stresses are evaluated assuming infinitesimal deflections and by use 


of the material characteristics of the finite elements. 

The stiffness finite element approach or approximation to the true struc- 
ture provides displacements which are a lower bound on the average. Lower 
bound on the average means that the true average deflection of the mirror is 
never less than the computed average deflection. The true solution is con- 
verged upon as the finite elements are reduced in size. 

In addition to the verification of the accuracy of ELAS, a verification of 
the adequacy of the finite element size is necessary since the results provide a 
nonconservative lower bound. In many situations the mirror was subdivided 
and the significant subdivisions were analyzed separately. This allowed more 
detail evaluation because of the ELAS limitation in the maximum allowable 


nodes or degrees of freedom. 

Check Problems 

As indicated in the previous paragraph, the adequacy of ELAS and mesh 
size was verified. 

First, an infinitely long thick cylinder with a thermal gradient through 
the thickness (see Figure 5) was solved with ELAS. The problem was chosen 
because an exact solution is available from the literature (Theory of Elasticity 
by Timoshenko and Goodier). As shown in Figure 6, the correlation between 
ELAS and Timoshenko was within accuracy of the graph. 

The second problem was a rocket nozzle, in which an exact solution does 
not exist. A cross-section of the nozzle showing the computer elements and 
nodes (mesh) can be seen in Figure 7. However, an approximate solution to 
this complex geometric shape with thermal loading had been solved by use of a 
well established computer program referred to as "Rohm-Hass". As shown in 
Figures 8 to 15, a good correlation was obtained. 
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TEMPERATURE, °F 


Figure 5 



Infinitely Long Cylinder with a' Thermal Gradient through the Thickness 





r, in. 


SAMPLE PROBLEM 


Figure 6. Curve Showing Correlation of "Elas" 


fSSroof frawss l 



Program 
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Figure 7. A Section of Rocketdyne Engine Nozzle 
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Figure 8. Radial Deflection of the inner Surface for Rocketdyne Engine Nozzle 
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Figure 9. Radial Deflection of the Outer Surface for Rocketdyne Engine Nozzle 
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Figure 10. Thermal Deflection of the Inner Surface in Y Direction 
for Rocketdyne Engine Nozzle 
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Figure 11. Thermal Deflection of the Outer Skin in Y Direction 
for Rocketdyne Engine Nozzle 
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Figure 12. Longitudinal Section of Rocketdyne Engine Nozzle 
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NOTE: Rocketdyne Program Computed the Stress on each Elemeiit; while JPL Program Computed 
Stress at Node Points, The JPL Ploting Point is the Average of Four Corners of each 
Element, 


Figure 13, Rocketdyne Engine Nozzle, Section A-A 
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Figure 15. Rocketdyne Engine Nozzle, Section C-C 
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The solutions of samples provided confidence in ELAS. 

EXAMPLE PROBLEM 

A typical example is the evaluation of the thermal distortion of a mirror 
subjected to a 3°F rise from inside to outside in the fluid used to cool the mir- 
ror. As shown on Figure 16, the fluid passed through the cooling holes in the 
center of the mirror. Because of geometric and loading axisymmetry, only 
one-half of any mirror cross-section adequately defines its structural behavior. 

As shown in Figure 16, the mirror is idealized by 247 nodes and 190 finite 
elements. The temperatures for each element were obtained from the temper- 
ature control computer program defined in Appendix C, 

A sample computer output describes in detail the necessary input and the 
available output. The explanation of significant values are: 

(1) Table 8 with the title of the analysis GREGORIAN PRIMARY 
MIRROR summarizes the control parameters. 

(2) Table 9 titled MATERIAL PROPERTIES , defines the material 
properties of the U.L.E, fused silica used for the analysis. Note 
that the properties are assumed constant over the temperature 
range for this particular problem. Also in Table 9, TEMPERA- 
TURE INCREASE TYPES defines the total different temperature 
increases in the finite elements used for the analysis. 

(3) Table 10 defines NODAL COORDINATES as the X and Y locations 
of each node. . 

(4) ' Table 11 defines the MESH TOPOLOGY as the definition of each 

finite element. Each finite element is completely defined by the 
node points which define the element as well as the material char- 
acteristics and temperature changes. 

(5) Table 12 gives the NODAL DEFLECTION output. For each node 
the displacement along the X and Y axes are given. 
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Table 8, Linear Elasticity Problem 


GREGORIAN PRIMARY M IRKOR 


TOTAL NUMBER OF NODES 247 

TOTAL NUMBER OF FINITE ELEMENTS 190 

DEGREES OF FREEDOM AT A NOCE 2 

ITYPE VALUE 0 

IGEM .VALUE 0 

ISTR 'VATUE 0 

MAXIMUM NUMBER OF CONTACTS IN AN ELEMENT 4 

CONTACT NUMBER INCLUDING DUMMIES , 7 

IBN VALUE 10 

TOTAL NUMBER OF CONCENTRATED LOADS 0 

PRESSURE TYPES 0 

MATERIAL TYPES I 

THICKNESS TYPES 0 

TFFFSraturE Chang'e 1 types .99 

TEMPERATURE GRADIENT TYPES ALONG Y ~0 

TEMPERATURE GRADIENT TYPES ALONG l -0 

AREA TYPES 1 -0 

TO KSTOTTO NTOfFT TYPES =V 

Y MOMENT OF INERTIA TYPES ”0 

ZHTOTTETirTOF INERT Ii tVpSs =0 

NUMBER OF ANGLES FIXING PRINCIPAL AXES -0 
I NX VALUE 4 

INP VALUE 1 

— r 

ICOR VALUE . -o 

— toottvatue ^ 

IMES VALUE -0 

— tpt k " . v &t urrmr~s'H err wcal nodal axes - -o' 

CHAIN PROGRAM TAPE NUMBER 2 

'SCRATCH' TAPE' NUWfTE W 

acceleration»unit MASS -0. 

— directicn cosines op 'A'C' c'ErERxrnTfi ^or 


0 FOR ISOTROPIC, 1 FOR ORTHOTROPIC, 2 FOR GENERAL MATERIAL 

0 FOR 2-, 1 FOR 3-DIMENSIONAL STRUCTURES 

1 FOR PLANE STRAIN CASE, 0 OTHERWISE 


NUMBER OF SUPRESSEO DEGREES OF FREEDOM IF NO MULTIPLE CONNECTIONS 


NUMBER OF THE LINK FROM WHICH RETURN-TO-BEGINNING IS MAOE 

0 MINIMUM PRINT, 1 PARTIAL PRINT, 2 COMPLETE PRINT 

0 NO REIAeELLING, 1 RELABEL, 2 OR 3 READ CARDS FOR RELABELLING 

0 READ CARDS, 1 CALL SUBROUTINE COftG FOR COORDINATES 

o read Cards, i CALV'TObroutine bung for boundary conditions 

0 READ CARDS, 1 CALL SUBROUTINE MESG FOR MESH TOPOLOGY 

0 ASSUME TEKTTrT COMPUTE AS (• R iNCTfOSTJ ~T «'£ ad as input 

0 OR - 00 NOT COMPOT'£ - R ESlDUALS, OTHERWISE COMPUTE 



Table 9. Material Properties 


MATERIAL PROPERTIES 
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Table 10. Nodal Coordinates 
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Table 10. Nodal Coordinates (Contd) 
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Table 11. Mesh Topology and Element Property Types (Contd) 
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0 

0 

0 

0 
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0 

1 

0 
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0 

“0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

16 
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1 
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-0 
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0 

0 

0 

0 
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0 

0 

0 

c 

16 

c 

1 

0 

93 

1 0 

-0 

0 

0 

0 

0 

0 


130 
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165 

0 

0 

0 

A 

V 

16 

0 

1 

0 

33 

0 

-0 

0 

0 

0 

0 

0 


131 

165 

172 

173 

166 

0 

0 

0 

c 

16 

0 

1 

0 

35 

0 

-0 

0 

0 

0 

0 

0 


132 

166 

173 

174 

16? 

0 

0 

0 

0 

16 

0 

1 

0 

48 

0 

-0 

0 

0 

0 

0 

0 


133 

167 

174 

175 

163 

0 

c 

0 

0 

16 

0 

1 

0 

68 

0 

-0 

0 

0 

0 

0 

0 


134 

168 

175 

176 

165 

0 

c 

0 

0 

16 

0 

1 

0 

83 

0 

-0 

0 

0 

0 

0 

0 


135 

lo9 

176 

177 

170 

0 

c 

0 

0 

16 

0 

1 

0 

95 

0 

-0 

0 

0 

0 

0 

0 


136 

171 

178 

179 

172 

0 

0 

0 

0 

16 

c 

L 

0 

34 

0 

-0 

0 

0 

0 

0 

0 


137 

172 

179 

180 

173 

0 

0 

0 ' 

G 

16 

0 

1 

0 

26 

0 

-0 

0 

0 

0 

0 

0 


138 

l 74 

18 i 

18 2 

1/5 

0 

0 

0 

G 

16 

0 

1 

0 

65 

0 

-0 

0 

0 

0 

0 

0 


139 

175 

18 2 

183 

176 

0 

0 

0 

0 

16 

0 

1 

0 

82 

0 

-0 

o 

0 

0 

0 

0 


14C 

176 

183 

184 

1 77 

0 

c 

0 

c 

16 

0 

1 

0 

95 

0 

-0 

0 

0 

0 

0 

0 


141 

rnr 

183 
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179 

0 

c 

0 

0 

16 

0 

1 

0 

39 

o 

-0 

0 

0 

0 

0 

0 


142 

179 

136 

137 

180 

0 

0 

0 

0 

16 

0 

1 

0 

41 

0 

-0 

0 

0 

0 

0 

0 


iv 5 

130 

187 

183 

181 

0 

0 

0 

0 

16 

0 

1 

0 

53 

o 

-0 

0 

0 

0 

0 

0 


144 

181 

18 8 

189 

182 

0 

0 

0 

0 

16 

0 

1 

G 

7 L 

0 

-0 

0 

0 

0 

0 

0 


145 

132 

189 

190 

183 

0 

0 

0 

0 

16 

0 

1 

0 

87 

0 

-0 

0 

0 

0 

0 

0 


146 

183 

190 

191 

184 

0 

0 

0 

0 

16 

0 

1 

0 

96 

0 

-0 

0 

0 

0 

0 

0 


"T4'7 ‘ 

1ST 

192 

L93 

186 

0 

0 

0 

0 

16 

0 

1 

0 

39 

0 

-0 

0 

0 

0 

0 

0 


148 

186 

193 

194 

187 

0 

0 

0 

c 

16 

0 

1 

0 

33 

0 

-0 

0 

0 

0 

0 

0 


149 

188 

L95 

196 

189 

0 

0 

0 

0 

16 

0 

1 

c 

69 

0 

-0 

0 . 

0 

0 

0 

0 


. 150 

189 

196 

197 

19C 

0 

0 

0 

0 

16 

0 

1 

0 

86 

0 

, -0 

0 

0 

0 

0 

0 


— rsr - ■ 

— FJtT 

1<T7” 

198 

191 

o 

0 

0 

c 

16 

0 

1 

0 

96 

0 

-0 

0 

0 

0 

0 

0 


152 

192 

199 

2C0 

193 

0 

0 

0 

0 

16 

0 

1 

0 

45 

0 

-0 

0 

0 

0 

0 

0 


V5~3 ~ ' 

1 T 3 ~ 

T 37 T - 

?0L 

194 

0 

0 

0 

0 

16 

0 

1 

0 

48 

0 

-0 

0 

0 

0 

0 

0 


154 

194 

201 

202 

19b 

0 

0 

0 

0 

16 

0 

1 

0 

57 

0 

-0 

0 

0 

0 

. 0 

0 


pj-5 

L95 

202 

203 

196 

0 

0 

0 

0 

16 

0 

1 

0 

76 

0 

-0 

0 

0 

0 

0 

0 


156 

L96 

?0i 

204 

197 

0 

c 

0 

0 

16 

0 

1 

0 

90 

o, 

-0 

0 

0 

0 

0 

0 


T5 - 7 

j£V? ' 

204 

T05 

198 

0 

0 

0 

c , 

16 

c 

1 

0 

97 

0 

-0 

0 . 

0 

0 

0 

0 


158 

199 

206 

207 

200 

0 

0 

0 

0 

16 

0 

1 

0 

44 

0 

-0 

0 

0 

0 

0 

0 




Table 11. Mesh Topology and Element Property Types (Contd) 


rs9 

~2W 

—i< rr 

roa 

201 

0 

0 

0 

0 . 


16 

c 

1 

0 

38 

'0 

-0 

0 

0 

0 

0 

0 


160 

202 

20 9 

210 

203 

0 

0 

0 

0 


16 

0 

1 

0 

72 

0 

-0 

0 

0 

0 

0 

0 


16 L 

~TST 

2i<r 

?ll 

204 

0 

0 

0 

0 


16 

0 

1 

0 

88 

0 

-0 

0 

0 

0 

0 

0 
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204 

211 

712 

205 

0 

0 

0 

0 


16 

0 

l 

0 

97 

0 

-0 

0 

0 

0 

0 

0 


“ - T&3 

? 0 6 

2L3 

214 

207 

0 

0 

0 

c 


. 16 

0 

1 

0 

So 

0 

-0 

0 

0 

0 

0 

0 


164 

207 

214 

215 

208 

0 

c 

0 

u 


16 

0 

1 

0 

52 

0 

-0 

0 

0 

0 

0 

0 


res 

2or 

213 

2T5 

209 

0 

0 

0 

0 


16 

0 

1 

0 

60 

0 

-0 

0 

0 

0 

0 

0 
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216 

217 

21 C 

0 

0 

0 

0 


16 

0 

1 

0 

80 

0 

-0 

0 

0 

0 

0 

0 
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2 10 

2H 
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2Ti 

0 

0 

c 

c 


16 

6 

L 

0 

93 

0 

-0 

0 

0 

0 

0 

0 
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0 

0 

0 

c 


16 

0 

i 
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0 
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0 

c 
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1 

0 
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0 


m 
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2D 

0 
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15 

c 
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0 

1 

Q 
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-0. 

0 

0 

0 
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0 
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0 

0 
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0 

1 
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-0 

0 

0 

0 

0 

0 


173 

218 

225 

2 26 

219 
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0 
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0 
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0 

0 


183 

ITT 

238 

~TV) 

T32' 
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0 
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-0 
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0 
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0 
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0 
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0 
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0 
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AT 0. 

SEC. 

OF RELABELLING 

-0 

SHEEPS 

PERFORMED . 

UPPER 

0FF-8AN0 ELEMENT 

COUNT OF 

MESH 

TOPOLOGY 

MATRIX 

IS 

27390 



AT 8.95 

SEC. 

OF 8 ELA8ELL ING 

370 

SHEEPS 

PERF0RM60* 

UPPER 

OFF-BAND ELEMENT 

COUNT OF 

MESH 

TOPOLOGY 

MATRIX 

IS 

28305 





Table 12. Nodal Deflections 


INPUT LINK TlCis 24.42 
GENERATION LINK 'f OCR 44 

SECONDS. 

.30 SECONDS. 







NODAL 

DEFLECTIONS 




NODE 

DISP. ALONG X 

CISP. ALONG Y 

DISP. ALONG Z 

ROTA. ABOUT X 

ROTA. ABOUT Y 

ROTA. ABOUT Z 

L 

0 . 

C.2092358F-06 

C. 

0. 

0. 

0, 

"2 

0 . 

o.-ysrrnos^cir 

0. 

0. 

0* 

0. 

3 

0. 

0.2758441E-C6 

0. 

a. 

0. 

0. 

4 

57 — — 

0.29L6618E-C6 

c. 

0. 

0. 

0. 

5 

0 . 

0. 304 1 5 79 E- 06 

0. 

0. 

0. 

0. 

6 

0. 

0. 3165377 E-C6 

G. 

0* 

0. 

0. 

7 

0.155 /CCOt-06 

0. 19540P3E-C6 

0. 

0. 

0. 

0. 

3 " 

0. 1S8C57SE-C6 

0.241 976P8-C6 

c. 

0, 

o* 

0. 

9 

0 . 1 738536£-06 

0 *266 70 33 £-”06 

0. 

0. 

0. 

0- 

10 

0.1770 32? 6-C 6 

0. 2017393 E- 06 

0. 

0. 

0. 

0. 

11 

0. l7Vfc670c-C6 

0.29559508-06 

0. 

0. 

0. 

0. 

12 

0 . 151 IC788-C6 

0.30d3?03E-06 

0. 

0. 

0. 

0. 

L3 

0. I996930E-96 

0.18574408-06 

0. 

0. 

0. 

0- 

14 

0 .2136473L-C6 

0.231b363E-06 

0. 

0- 

0. 

0. 

15 

0 .2275630E-0O 

0.25747558-06 

0. 

0. 

0. 

0. 

16 

3TT32 29 54f7cfc 

0.27270 lOE-CL? 

0. 

0. 

0. 

0. 

17 

0 .236J394t-C6 

0.2074685E-C6 

0. 

0. 

0. 

0. 

13"*” 

0 .23999C2E-t6 

0.301 1412E-C6 

0- 

0. 

0. 

0. 

19 

0. 

-9.45145678-Of 

0. 

0. 

0. 

0. 

20 

0. 

0. 3213373E-08 

0. 

0. 

0. 

0, 

21 

0.19 10 199E-G6 

0.10154508-06 

0. 

0. 

0. 

0. 

22 

0 .2215747E-0O 

0. 1791909 8- C6 

0* 

0. 

0. 

0. 

23 

0 .24473728-06 

9. 22455368-06 

0. 

0. 

0. 

0, 

24 

0 .25483208-06 

0.24954 858-C6 

0. 

0. 

0. 

0. 

25 

0.2606476b-06 

0.2659079 8-C6 

0. 

0. 

0. 

0. 

26 

0 .266888 3C— 06 

0.28103438-06 

0. 

07 

0. 

0. 

27 

0.27059148-06 

0 . 2969502 E-06 

0. 

0. 

0. 

0. 

28 

0. 

0. 11734958-06 

0. 

0. 

0. 

0. 

29 

0 . 146o822E-06 

-0.87223008-07 

0. 

0. 

0. 

0« 

30 

0.2 1982 17E-06 

9.037051 68-C 7 

0. 

0. 

0. 

0, 

31 

0.2643129 £-06 

0.1593299E-06 

0. 

0. 

0. 

0. 

32 

0.297634SE-06 

0.20375738-06 

c. 

o. • 

0. 

0. 

• 33 

0.3111C43E-06 

0.22870638-06 

0. 

0,. 

0. 

0. 

34 

0T3'f8r6'£7E~C6" 

" *" 0 . 2*45 52T2C-C 6 

c. 

0. 

0. 

0. 

35 

0.32487576-C6 

0.2636970E-C6 

c • 

0. 

0. 

0. 

36 

0. 1212834S-C6 

0.28130C4E-C6 

0* 

0. 

0. 

0. 

37 

0.14405218-06 

-0.667C2 70E-09 

0. 

0. 

0. 

0. 

38 

OTIJ'SOCo/E-'JB 

0. 1086891 

0. 

0. 

0. 

0. 

39 

0 . 24 3 3 460E-C6 

9- 59409 1 6E— 07 

0. 

0. 

0. 

0. 

40 

0.39*5400Lt-'06 

9.12925468-06 

0. 

0, 

0. 

0. 

41 

0. J524555E-06 

0.17171428-06 

0. 

0. 

0. 

0. 

42 

0 .370 5931 t-06 

0. 1957322 06 

0. 

0* 

0. 

0. 

43 

0. JB43449E-06 

0.2 13049 3E— 06 

0. 

0. 

0. 

0. - 

44 

<5T3v6n04 Lc-06 

0 6 

0. 

0. 

0* 

0. 

45 

0 .40263348-06 

0.24566648-06 

c* 

0. 

0. 

0. 

4 6 

0 .18.494 r2*H77(T5 _ 

— =95T8'5T6T9FE=r7" 

c. 

0- 

9. 

0. 

47 

0.10561.498-06 

0. 

Q . 

0* 

0. 

0. 

48 

0.26576498-06 

0.4690699E— C7 

0- 

0. 

0. 

0. 

49 

0 .326 1? 538-06 

0.10879446-06 

c. 

0. 

0. 

0. 

5o 

0.380780Yfc-<5*6 

0.15018328-06 

G. 

0. 

0. 

0. 

51 

U .40533748-06 

0.17242438-06 

G. 

0. 

* rj 

0. 

0. 



Table 12. Nodal Deflections (Contd) 


52 

0*418 i 8906-06 

0.10962302-06 

0 . 

0 . 

0 . 

0 . 


53 

0.6 32 20696—06 

0.21018 17E-06 

0 . 

0 . 

0 . 

0 . 


54 

0.66 5 1 1 32S'-Vj6 _ 

0 • 23201 G 1 E—06 

0 . 

0 . 

0 - 

0 . 


55 

0 .17338352-06 

0.10626292-06 

0 . 

0 * 

0 * 

0 . 


?6 

0*<i 1009 57 E-04 

-Ol'B'693999-f-CT' 

0 . 

0 . 

0 . 

0 * 


57 

0.20731 168-06 

- 0.1253610 2- C7 

0 . 

0 . 

0 . 

, 0 . 


ra 

7r.T06“)TS3 2^5 1 

iT.lSMO 3 82- (77 

0 * 

0 . 

0 . 

0 . 


59 

0.297 LO44E-06 

— 0 • 1265 142E-07 

0 . 

0 . 

0 * 

0 . 


60 

C.377JJ62F-06 

0 . 6 5590072-C7 

"o* 

0 * 

0 . 

0 . • 


61 

0.669112 9 £-06 

0 *8716688 E-C7 

0 . 

0 . 

0 . 

0 . 


62 

0 .67901781-06 

0.10B63672-C6 

0 . 

0 . 

0 * 

0 . 


63 

C. 699 35506-06 

0.12586522-06 

0 . 

0 . 

0 . 

0 . 


64 

0.517o002f-C6 

0.16656692-06 

0 - 

0 . 

0 . 

0 - 


65 

0 .53 372 786-06 

0. 16963682-06 

0 . 

0 . 

0 . 

0 . 


66 

0*3125 19 7F-06 

-0. 60821202-07 

0 . 

0 - 

G. 

0 . 


67 

0.39669352-06 

0. 13833562-C7 

0 . 

0 . 

0 . 

0 . 


68 

0.678 11 1 12-06 

0.55970C7E— 07 

c. 

0 . 

0 . 

0 . 


69 

0. 51329732-06 

0.77271352-07 

0 . 

0 . 

0 . 

0 . 


70 

0*5336811 E- 06 

0.9S282296-C7 

0 . 

0 . 

0 . 

0 . 


71 

0.55370652-06 

0. 1 1617502-06 

0 . 

0 . 

0* 

0 . 


72 

0.37182652-06 

0.13839692-06 

0 . 

0 . 

0. 

0. 


73 

0 * 352099 lt -06 

-0.12677562-06 

3. 

0. 

0 . 

0 * 


74 

0.66699716— C6 

-0 . 7 56 88 6 12-07 

0 . 

0 . 

0 . 

0. 


75 

0. 56677652-06 

-0.32622832-07 

0 . 

0 . 

0 . 

0. 


74 

0 .5889390^7^6 

-5. lOTSTlTE^O 7 

0 * 

0 . 

0 * 

0 . 


77 

0 .0 16676 3 2-06 

3.82175P 12-08 

0 . 

0. 

0 . 

0 * 


78 

0 *64069 3*1 E-C 6 

0.29301822-07 

G. 

0. 

0. 

- 0 . 


79 

0.66161122-06 

0.51723212-07 

0 . 

0 . 

0. 

0. 


80 

0 .36059732-06 

-0. 165092 12-06 

0 . 

0 . 

0* 

0 . 


81 

0.47071 116-06 

-9-1 166035 E— C 6 

0 . 

0 . 

0 . 

0 . 


82 

0TT75m7692-06 

-0.7 36907 7E-07 

0 . 

0* 

0* 

0. 


83 

0 .6236 357 2-06 

-0 • 5 1 56680E-C 7 

0 . 

0 . 

0 . 

0 . 


86 

0.66969272-06 

-0.33866552-07 

0 . 

0. 

0. 

0. 
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APPENDIX B 

EXPLANATION OF ONE-HALF OPTICAL PATH 
DIFFERENCE CALCULATION 

To calculate the l/2 O.P.D, the following analysis was used. First, a 
mathematical parabola of 250 cm focal length was derived. The equation for 
this is 

y » X 2 /4f = 0.00254000002X 2 in.; f = 98.425196 in. 

Then the computer distortion data Ay was listed for selected X radii. Then the 
error in distortion (Ay^) at X = 0 was forced to zero by making y = 0 at 
X = 0. Then each corrected error Ay - Ay^ was computed. The new surface 
of revolution was described mathematically by letting y" = y + (Ay - Ay^). To 
determine how closely the new curve, y", approximated a new parabola (thus a 
new optical surface of satisfactory imaging characteristics), a new mathemat- 
ical parabola (y) was generated, choosing a focal length (f ) based on the 
y + (Ay - Ay^) value found at radius X of 12,79 inches (32.5 cm) from the center 
of. the mirror. (See sample calculation in Table 4.) It was found through 
experience that this focal length gave the least error (difference along the y 
axis) between the new parabola (y) and the y" curve. The sum of the maximum 
positive and negative deviation between the two curves (y'« y") will give the 
l/2 O.P.D. In the sample calculation in Table 4 this amounted to 1. 05 micro- 
inches or a l/2 O.P.D. of \/l9. 
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APPENDIX C 

THERMAL ANALYSIS OF THE PHOTOHELIOGRAPH PRIMARY MIRROR 

The initial primary mirror configuration (a uniform thick disk) was 
established for the thermal investigation. The diameter of the disk was estab- 
lished by the current dimension of the mirror at that time which was 80 cm 
(31. 5 in. ), with a thickness of one inch (see Figures 17 and 18). 

The boundary conditions for the model were: 

1. The front surface was sun illuminated with the exception of a shaded 
circle in the center which had a diameter of 10 cm (3. 94 in. ). 

2. The back surface had a uniform temperature of 70 °F. 

3. The edges were adiabatic. 

As usual, the model was discretized into convenient elements, for the 
purpose of the computer analysis. Figure 17 shows the model which was 
programmed. 

Once the mathematical model was completely described, a program 
could be written to do the required analysis. The analysis was accomplished by 
use of finite difference techniques. These techniques were implemented by 
means of a general heat transfer computer program named Chrysler Improved 
Numerical Differencing Analyzer, commonly known by its acronym, CINDA. 

This program has the capability to do the thermal analysis as well as the 
menial calculations necessary to set up the thermal network. By listing mate- 
rial properties', and geometrical characteristics a subroutine can be called on 
to multiply the volume of a node by the density and specific heat and arrive at 
the thermal capacity of the node. In a similar manner, constant values can be 
loaded and called up in a subroutine to calculate the thermal conductance 
between a pair of nodes. 

To illustrate the program, a sample problem is presented and explained. 
In Table 13, the first card titled BCD 3THERMAL SPCS prepares the com- 
puter system for a thermal problem. The following three cards are title cards 
describing the boundary conditions. 

The network is set up by the card "BCD 3NODE DATA" and the cards that 
follow down to the card "END". Note that this data is listed in triplets. The 
first data in the triplet lists the node numbers, the second data value is the 
initial temperature of the node, and the third data value is the capacity of the 
node. In the particular example shown, the capacitance is zero. Later in the 
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NOTE DISTORTION TO SCALES: 

RADII TO CENTERS OF TOP ROW OF NODES 

r 70 = ° r 80 = h0 In * r 90 ~ 1,93 in * r 2 = 2 * 75in * 

r = 3.93 in. r. = 5.50 in. r_ = 7.07 in. r = 8.64 in. 

3* 4 o o 

r 7 = 1 0.2 in. = 1 1 .80 in. r ? = 13.35 in. r 1Q =14.50in. 

^ = 15.30 in. 

DEPTHS TO CENTERS OF FIRST COLUMN OF NODES 


7 q = 0. 05 in. 

X 71 = 0.20 in. 

X^ = 0.40 in. 

X ?3 = 0.60 in. 

74 = 0.08 in. 

X-... = 0.95 in. 
/o 




Figure 17. Subdivision of Plane Gregorian Mirror for Heat Transfer Calculations 



TEMPERATURE 



RADIAL POSITION, In. 


Figure 18. Steady State Temperature Distribution in Gregorian Mirror 
with 10cm. dia. Shaded Spot and Uniform 70°F Back 
Surface Temperature 
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Table 13. Chry.sler Improved Nume 


i 

i 


9DVEKALL ANALYSIS OF TAPERED ORE COR I AN MIRROR rfITH 1NTE 

- *RH"A'rx’o3'L iwsr ctwnme cth at l-u L£'“a o ar ti 

5FL0H -TURB., UNIFORM IN TO OUT 

3N00E DATA 

TTTOTrff7ZT7&~0T3 C7 DC. 37*77 OV 7075 (TO DCS C/O'. ,‘0“ 


7,70. , 0, B , /0. ,0,9 

■ I SmJT,7T,V477 o . , o' 

19. 70. . 0.20. 70. »C 

25, /o:, 0,26, 70., 0 

31. 70. . 0.32. 70. .0 
J777C7T0V3 S ; 7 O'C, 0" 
43j70. ,0,44,70.,0 

— 49 ,"70'. ,'0','5 O', 7'0 . , 0 . 
55^70. ,0,56,70. ,0 
"61 , 70". (0 ,'62, 70."* 0 

67. 70. . 0.6d,70. ,0 
— 73T7 OV," 0774 ,'7 0 .' CO 

79.70. . 0.80.70. .0 
- *'8b", 70. , 0( 86, 70. , O' 

91.70. . 0.92.70. .0, 
' 97, 70.', 0,98, 70. ,0. 


103.70. , 

— r&vrror, 

115.70. , 
~121 ,70.', 

127.70. , 

I3'j,7'0."i 

139.70. , 
17377(37, 

151.70. , 
“T57, 70." | 

163.70. , 
76*9, 7 OV, 

175.70. , 
T6T,T(n"i 
187,70-j 

' 191 ,70., 
_iy5,70. , 
199, 70 (( 
203, 7 0_._, 
2(77770. , 
_211,70._i 
"215,"7*0. , 

219.70. , 
' 2 23 ,~7 0 . , 

227.70. , 
“Z5T77B71 

235.70. , 
*239 ( 70."', 

243. 70. 
■'-2'47, 70. 

-252,70. 

- 257770 . 


, 70., 0, 10, 70. ,3, 11,70., 0, 12,73.,0 
,'l“5770'.“,0 ,16', 73".',0,i7 , 7 7 . , O', 18'," 7 0.(0 
,21,70. ,0,22,73. ,0,23,73. ,0,24, 70., 0 
,27,70. ,0,28", 73. ,0,29,70. ,0,30, 70., 0 ' 
,33,70. ,0,34, 73. ,0,35,70. ,0,36,70., 0 
,39,73. ,0, TO , 73 . , 0 , 41 , 73., 0,42, 70. ,0 
,45,7 0 . ,0,46,70. ,0, 47 , 70., 0,48, 70., 0 
,'51,70. ,0,52,7 0. ,0,53,70. ,0,54, 70., 0 
,57, 70. , 0,58, 70., 0,59, 7 3., 0,60, 70., 0 
,63,70. ,'0, 64, 70 (, 0,65 ',7 3. ,0, 66, 70., 0 
,69,70. ,0,70,70. ,0, 71,70. ,0,72, 70., 0 
,'7 5', 7 OV , O'T T'67737T0 , 77’, 7 3 .' ,'0", 78', 70.', 0 
,81,70. ,0, 82, 70., 0,83, 70., 0, 84, 70., 0 
,87,70. ,0,88,70. ,0,89, 7 0., 0, 90, 70., 0"““ 
,93,70. ,0,94,70. ,0,95,7 3. ,0,96, 70., 0 
,99, 70. ,0,100(70. ,0,101, 70. ,0, 102,70. ,"(T 


0 , 104 ,70 . , 0 
"0(110", 70. ‘,0 
0, 116, 70. ,0 
0, 122 ,*7C * , 0 
0, 128, 70. ,0 
0,134, 70., 0 
0, 140j.ro., 0 
"(T, 146,70. To 
0,152,70.,0 
0( 158,70. ,'0 
0,164,70. ,0 
0,170, 70.", O' 
a, 176,70. ,o 
‘oTleZVVo. , q 
0, 188 , 70. , 0 
-TT, 1*92 ,"70." 
-1., 196, 70. 

-T; ,'2-00, 70 : 

- 1 . ,20 4, 70_. 
-1 . ,208(21). 
“ 1 ., 212 , 73 . 
‘-£., 216 , 70 . 
-1., 220, 70. 
- 1 . , 224 , 73 ( 
-1., 228, 70. 

-rrr232,rcr." 

-1., 236, 70. 
-1.(240, "70. 
— 1 . , 2 4 4_, 70. 
( 0 ( -2 4*8 »*7" 3 » 
,0,-253,70. 
7D,-2"5"a7?ov 


LOS, 70., 0,106, 70. ,0,107 
L'1'1 ,"70 . , 0,"1 1 2V"70". ,0,113' 

117. 70. . 0.118. 70. .0.119 
L23, 70. ,0,124, 70. ,0,125 
L29, 70. ,0,130, 73. ,0,131 
135 , 70. ,0,136, 7 3 VC 0 ," 137 
14L,70.,0,1_4 7, 70^,0,143 
147 ,~76'. , O', 1 4l>" ,70'. ,"0 ,"l 49 

153. 70. . 0.154. 70. .0.155 
1 5 9 ,T OV, 0 ( 1 S'OV'7 0'.70‘( 161 
16?, 70., 0,166, 70. ,0,167 
171, 70 V, O', 1 7 2 , 70V, 0*, f 7 3 

177. 70. . 0.178. 70 0, 179 

’ ra'3','70 . , 0 , f3~, 7 O'. , 0 , 1 B 5' 

189. 70. . 0.190. 70. .0_ 
-l'C, 19*3, 70. , -1. ,’£94770. 
-1., 197, 70*, -1., 19 8, 70. 
-1. ,20T,70".,-1.,"20Z", 70. 
-1. ,205,70. , -1., 206, 73. 
-TT72"0 9T'7" STT^l” , 2"i o"77 o . 
-i.,213,70.,-1.,214,70. 
- 1( , 2"1 7","70 . ( - 1 . , 21 8 ,'70 . 
-1. ,221, 70., -1., 222, 70. 
”1. ,225, 70. , — 1* ,226, 70. 
-l.,229,73.,-l.,230,70. 
-1. , 233T7"0. , -l"( , 23-4T/0: 
-1. ,237, 70.,-!., 238, 70. 
-1. ,241,70.,— 1*, 2*42 , 70. 
-1. , 245, 70. ,-l., 246,7'). 
0(-"249, 70. ,0,-260," 73. ,0 
0,-2 54 , 70. ,0,-2 5 5 , 70 . , 0 
"0.-2T977O. , 3,-250, 70. , 0 


70.. 0. 108. 70. 
"70., 0, 114,70. 

70. . 0. 120.70. 
70".C0,126( 70". 

70.. 0. 132. 70. 

70. . 0. 138. 70. 

70. . 0. 144. 70. 

7 o'.'Co CTs'o'Cvo". 

70. . 0. 15., 70. 
'70".,0,"16""2,l0. 

70. . 0. 168. 70. 
'70'., 0, 174", 70( 
7 0 ,_,_0 , 18 0, 70_. 
TO"., 0","l86"7TO. 

-lV 

- 1 . 

-'1.' 

_ 

_i._ 

-1. 

-1.' " 

-1. 

'- 1 . 

-1. 

~r. 

-1. 

-i.‘ 

-1. 

-251,7u.(0" 
-256, 70. ,0 


END 


al Differencing Analyzer 




Table 13. Chrysler Improved Numerical Differencing Analyzer (Contd) 


RELA TIVE ACTUA L 

NODE NUMBERS “ 'JUDE NUMBERS 


1 

THRU 

10 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


11 

THRU 

20 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


21 

THRU 

30 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 


31 

THRU 

40 

31 

32 

33 

34 

35 

36 

37 

38 

39 

<*0 


41 

THRU 

50 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


51 

THRU 

60 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


61 

THRU 

70 

61 

62 

63 

64 

55 

66 

67 

68 

69 

70 


71 

THRU 

80 

71 

72 

73 

74 

75 

76 

77 

78 

79 

BO 


81 

THRU 

90 

. 81 

82 

83 

84 

65 

86 

87 

88 

89 

90 


91 

THRU 

100 

91 

92 

93 

94 

95 

96 

9 7 

98 

99 

100 


101 

THRU 

no 

101 

1 02 

103 

104 

105 

106 

10 7 

108 

109 

110 


11L 

THRU 

120 

111 

112 

113 

1 14 

115 

116 

117 

1 18 

119 

120 


121 

THRU 

130 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 


131 

THRU 

140 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 


141 

THRU 

150 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 


151 

THRU 

160 

151 

152 

L53 

154 

155 

156 

15 7 

158 

159 

160 


161 

THRU 

170 

161 

162 

163 

164 

165 

166 

167 

L68 

169 

170 


171 

THRU 

180 

171 

172 

173 

174 

176 

176 

177 

178 

179 

180 


181 

THRU 

190 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 


191 

THRU 

200 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 


201 

THRU 

210 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 


211 

THRU 

220 

211 

212 

213 

214 

215 

216 

21 f 

218 

219 

220 


221 

THRU 

230 

221 

222 

223 

224 

225 

226 

227 

228 

229 

230 


231 

THRU 

240 

231 

232 

233 

234 

235 

236 

?3 7 

238 

239 

240 


241 

THRU 

250 

241 

242 

243 

244 

245 

246 

24 7 

248 

249 

250 


251 

THRU 

260 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 



BCD 3CONDUCTOR DATA 

111.1 ,2,34,35,67,68,0 


112. 2, 3, 35, 36, 68, 59 , 0 

113.3, A, 36,37,69,70,0 
114, A, 6, 37, 38, 70,71,0 

115.5.6.38.39.71.72.0 

116.6.7.39.40.72.73.0 

117.7.8.40.41.73.74.0 

118.8.9.41.42.74.75.0 

119.9.10.42.43.75.76.0 
1110, 10, 11, 43,44,76,77", "6" 

1111.11.12.44.45.77.78.0 

1112', 'IT,' 1'3 ,45 ,46 V 7e',79To 

1113. 13.14.46.47.79.80.0 

1114. 14.15.47.48.80.81. 0 
1115, 15, L6,48, 49, 81,82,0 
1116, 16,17,49,50,82,83,0" 
1117, 17, 18, 50, 51, B3, 84,0 

1118. 18.19.51.52.84.85.0 

1119.19.20.52.53.85.86.0 
1120, 20, 21, 53, 54, 86, '87 , 0 

1121.21.22.54.55.87.88.0 

1122. 22.23.55. 56.86.69.0 
1123, 23, 24, 56, 57, B9, 90,0 

1124.24.25.57.58.90.91.0 

1125.25.26.58.59.91.92.0 
1 1 26726TT7 ,59,60,92,93,0 

1127.27.28.60.61.93.94.0 



Table 13. Chrysler Improved Numerical Differencing Analyzer (Contd) 


fTTBTTSTTO , 5X75279 r * . 9'5T0 
1129, 29, 30, 62, 63_, 95, 96,0 
' 'XT3'0T'3*OV3T763Tft,9'679>70 

1131.31.32.64.65.97.98.0 

l l'327727S'37ff57S67V8 , 99',£T 
L21, 1 ,39,39,67,0 

L'2T7273'573'S766,0~ 

123.3.36.36.69.0 

1-2-977737 r3-777 0 70' 

125.5.38.36.71.0 
r2"ff7S7"3'9T3'97 7"2" , 0 

127.7.90.90.73.0 

rr87-379n7rr,79-r0 

129.9.92.92.75.0 

rUlTJ , 1 0T9"37 9"3 i 767 0" * “ —7' 

1211.11.99.99.77.0 

72X27 1779 S" , 9 577870 

1213. 13.96.96.79.0 

12X4774 , 5T74TTfftJ7D ' 

1215.15.98.98.81 .0 

12X67X6, '99, W762,0 ‘ ~ 

1217.17.50.50.83.0 

rzT87TS;-5T73'i78T,'0' - 

1219. 19.52.52.85.0 

rZTOT^O", -5375-3786- 70 

1221.21.59.59.87.0 

1 2-?"27 2 27 5 5-755' , 88 , 0 " 

1223.23.56.56.89.0 
“ " 1229,29,57,57,90,0 

1225,25,58,58,91 ,0 

— nrzisT 2675 - 975-9 , 92 - 70 - 

1227.27.60.60.93.0 

im'tf 728-761761-, 99,0 “ 

1229.29.62.62.95.0 
~ 72307301 6 3 ,63 ,96,0 ’ 

1231.31.69.69.97.0 

I-2-32-7327S-576-579B-70' 

1233.33.66.0 

2233766. 99.0 

321.67.100.0 

‘ "-322768’, 101, O' 

323.69.102.0 

3T97T07 ro'37"0 

325.71.109.0 

326,72,1-91,0 

327.73.105.0 

32-8 774 7195,0 — 

329.75.106.0 

327077-6719976 

3211.77.107.0 

3212778,-20-3,0 

3213.79.108.0 

-"32T4, 8"0,"2-07 1 0' 

3215.81.109.0 
32 767827211,0 

3217.83.110.0 



Table 13. Chrysler Improved Numerical Differencing Analyzer (Contd) 


3^8,84, 21 3,0 

3219.85.111.0 

3220.86.219.0 " 

3221.87.112.0 
3222,88,223,0“ " 

3223.89.113.0 

3224.90.227.0 

3225.91.114.0 

3226.92.231.0 ~ 

3227.93. 115.0 
3228,94,23570 — 

3 229,95,116,0 

3230.96.239.0 

3231.97.117.0 

3232.98.243.0 

3233.99.118.0 
411 , 100,101, 0“ 

412.101.102.0 

413.102.103.0 
41 4, 103,104,0 

41~5, 104,194,6 

416.192.105.0 

417.105.198.0 
418,196,106,0 

419.106.202.0 

I 4110,200,107,0 

m 4111,107,206,0 

4112,204,108,0 

4113.108.210.0 

• 4114,208,109,0 

4115. 109.214.0 

4116.212.110.0 

4117. 110.218.0 
4118,216,111,0 

4119.111.222.0 

4120.220.112.0 
4121, 112, 226, "0 

4122,224,113,0 

4123.113.230.0 
4124,228,114,0 

4125.114.234.0 
4126,232,115,0 

4127.115.238.0 
4128,236,116,0 

4129.116.242.0 
4130,240,117,0 

4131, 117,246,6" 
4132,244,118,0 

421.100.119.0 
422,101, 1 20,0 

. 423, 102, 1*21, 6 

. ’ 424,103,122,0 

425.104.123.0 
426,193,124,0 

427.105.123.0 

428.197.125.0 




4, 12 . 5,0 
5X7,125, 126,0 


. 26, 

519,127,120,0 


O, 128,129,0 

5111,129,130,0 


5112.130.131.0 

5113.131.132.0 


LiirYV kktftl 


5115.1 33.13 4.0 
5116, 134, 135,0 

5117.135.136.0 


5127.145.146.0 
5128, r467r47T0" 

5129.147.148.0 
T1T571 48714 97CT 

5131.149.150.0 












Table 13. Chrysler Improved Numerical Differencing Analyzer (Contd) 


Tr577T5‘07r51T0' 

52_l,ll_9,U)3,0_ 

T22,120;i8^i0 

523.121.185.0 

■5277T227r5‘270 

525.123.153.0 

52'67T2T7rr470 

527.125.156.0 
*5787l’2*67 1*57 70 

529.127.158.0 _ 

3210,178,1*59,0 

5211.129.160.0 

T2T27XT07r5I75 

5213.131.162.0 • 

T2 177732 , 173 , 0 

5215.133.164.0 

T2T6Tf3'47Y6570 

5217.135.166.0 

■^TffTfTTTrsTTO 

_5_2 19,137 ,168 ,_0 

52*20, 138, 169**,0 

5221, 139,170,0 

5222.140.171.0 

5223.141.172.0 
“522777727X7370 

5225. 143.174.0 

"5 226, 744,175,0 

5227.145.176.0 
32187746, 7771*0 

5229.147.178.0 

5230.148.179.0 

5231.149.180.0 
*5*2*32 , 1 5*0, 1 81 , 0* 

3233.151.182.0 

611.183.184.0 

612.184.185.0 
"61*3,'l"85';i5'2,'0 * 

614.152.153.0 

615.153.154.0 

616.154.156.0 

617.156.157.0 

618.157.158.0 

6110. 158.159.0 

6111.159.160.0 

6112.160.161.0 

6113.161.162.0 

6114. 162.163.0 

6115.163.164.0 

6116.164.165.0 

6117.165.166.0 

6118.166.167.0 

6119 . 167 . 168.0 

6120.168.169.0 

6121.169.170.0 

6122.170.171.0 



Table 13. Chrysler Improved Numerical Differencing Analyzer (Contd) 


6124.172.173.0 

6125.173.174.0 

6126.174.175.0 

6127, 1 15,176,0 

6128.176.1 77.0 

6129.177.178.0 

6130.178.179.0 

6131.179.180.0 

6132.180.181.0 

6133. 181.182.0 : 

621.183.186.0 

622.184.187.0 _ 

623.165.188.0 

711,186, 187,0 _ _______ 

712. 187.188.0 

721.186.189.0 

722.187.190.0 

723.1 8 8.155.0 

811.189.190.0 

812.190.155.0 

901, 191 ,"247, 1*92, 24 7 » 193 , 247 , 194, 247, 0 

902.195.248. 196.248.197.248.198.248.0 
9Ui, l9V,24v,2WT249,201 ,249,20 2,249 , 0 

904.203.250.204.250.205.250.206.250.0 
9'0'5”',"2bT','25 1 , 2*08 , 25"f, 209 , 25*1 , 21 0, 251", 0 

906.211 .252.212.2 5 2.213 .2 52.214.25 2.0 

907.215.253.216.2 53,217,253,218, 25*3,0 

908.219.254.220.254.221.254.222.254.0 

909.2 23,25*57254,2 5*5,225,25*5,226, 255,0 
910, 227, 256, 228, 256, 229, 25b, 23D, 255,0 
911, 231, 23*7, 232,2*5 7*7233 ,"257, 23**4, 257,0 

912.235.258.236.258.237 .258.238.259.0 
913, 239, 259, 240*,'2*S9, 241, 25*9, 242/2*59, 0 
914,243,267), 244, 250, 245 ,260,246,260,0 

rfiii 

6CD 3C0NSTANfS DATA 

NLUCJP , rswri , 7063*5 , 27777 97*737 • 1*3*3 m*K,7.2*RH0,K3 = C 5'i)*B P 

4. 70. . 5. 0.43 6K4=INITIAL TtMP. , K5= ABSORBED SOLAR FLUX 
67T70. ,7,1 44 .* iK6=C0NVECTiVE COEFFICIENT 

10.0. , 12, 9 l, 13, 32, 14,0. ,.15,31,15,14,17,246,18,29 

Q— -jTj 

orlxca,o.ooi,arlx:a,o.ooi iTEMP. changes fur STEADY solution 
[97',*T3*0*,T9 87-*0O"57rt*97L7l*SK'l 9T= rfM £N D';KT93~ffU"T P UT7lU 9 9=*C S G'F'aT 

1301. . 800. 1302.. 523. 1303.. 310. 1304.. 719. 1305. .895 
1306, • 5*2*7 , f3*07 , 1. 43, 13*08,. 661,1 3*09, TTfa,~lfiO , 779*6"** 

13 11. 2. 08. 1312.. 931. 1313.2.41.1314.1.07. 1315.2.73 
l*3T67TT?07 nTT,*T.7JS7T3TSTI*r34Vr 3T973".*3'37 1 3 2 07 7747 

1321,3.70,1322,1.61 ,1323,4.03,1324,1.74,1325,4.35 

'1*3*2 6 V 1'. )f 771'3'2'774'r 5V7ri 2*87 2‘."'0 1 ,1*32975". 00,f3 30V 2-14 ** 

1331.5.33. 1332.2.28.1333.2.85.1400.1.36.1401. . 884 

1402.. 5*2 67* 1*4 0 3 ,* 17*2 2 ,“1*4*0 4* , 1*. 5 2 , 1 4 0 5V 2*. 4 2 , 1 4* 0 6 2 . 9 7 

1407,3.52, 1408,4.07,1409,4.62,1410,5.16,1411,5.72 
1412,6 .*2777713 ,*6 .*8*2 ,14L*47 7.3 7*, 14*15 , 77 9*0 V 1416*,8*.*47 '"** 

141 7, 9. 02, 1418, 4. 82, 1419, 2. 35, 1420, 1.75, 1421, 1.0 7 
14*22727*5?; 1423, 37 2*5 ,*1424, 1.9 5, 1425, 5. 46, 1*426, 2. 58 
1427, 7.03, 1428,3.22,1429,8.19,1430,3.57, 1431,8.75 



Table 13 . Chrysler Improved Numerical Differencing Analyzer (Contd) 


1*4327 3 . *8T,'14*3*379*756 , *143 4, 4*.*0 7 , 1*4 35 .'1 0 . 0 , 1 436*. 4.25 

1437, 10.5,1438,4.34,1439,10.5,1440,4.38, 1441, ID. 6 
r4V2 , 4*73*87*144 3,10. 5 , f4*4’4",4"73"lTl 4'45“, 10 “3, 1446 , 4 -"20 ' ” 

1447.9.89. 1448.4.01.1449.9.35.1450.3.80. 1451.4.49 

nT5T?74V4'07 i'4537S7L'37r454*, 278*5 *,X4*55 ; 27 0 3', 1'4 56 , 1 . 66 ~~ 

145 7,3.10, 14bH,7.8L, 1459, 3. 22, 1400,8.16,1461,3.56 
1^27Si r927r4'63T37r37l4'6'4'r9V49‘, 1V6574.0 4 ,*1 4 66 , 9794 ■■ 
1467,4.21, 1468,10.3,1469,4.29,1470,10.4,1471,4.32 
T7iTZ7T0747T4T374737)7r4T47Ti7n7r475',7.'Z 0',* 1*4*7 6*. 10. if 

1477.4.08. 1478.9. 55. 1479. 3. 86. 1480. 3.8 3. 1481 1 4. 49 
148*272 .'647178 3 , 570577484 ,*3730 i 1“ 4*85 , i:9 6, 14&6i S*.*24 
1487,3.42,1488,2.03,1489,5.25,1490,3.42 

rriT.~7i'S’,112,4.90,Yl3,5.'45',ll4,'4.7B,115W.20 

116,6. 06, 117, 8. 58, 118,8. 37, 119, 9. 29, 1110, 9. 76 

— rrrrrnTZTm 27irr4'm*ixrr2T57n'i*4Ti373*, 1 ns ,* 14 :i * 

1116. 15. 0. 1117.16.2.1118.16.2.1119.17.6.1120. 17.8 

: — ri'21 ;T9V57ii 22 ,797471 123 7217*2 ,'1124 ,2 1.0,1 125,22.8 

1126,22.7,1127,24.4,1128,24.3,1129,26. 1,1130, 26.6 

' '1X3172*77 8* ,'113*2 ,'4*3*. 5 

121.47.4.122.30.9.123.18.4.124.42.5.125.53.0 
l*2*573T.*T7T277ff4177r2Vr3 97T71 2*971 9V. 1'2*10T47*.*1” * 

1211.123.. 1212.55. 1. 12 13.142.. 12 14. 63. 1.1215.162. 
i'21'6'; 7*1*: 0 , 12 1 7 , 1 8 1 *12 L 8 ,*7*9*; 0, 12 19 ,‘200 . , 1*220 , 8 7 . 0 

1221. 219. . 1222. 95. 0.1223. 238. .1224. 103. .1225. 258. 

r22*SVm.'» *1*227,277*., ' 1220,T19.', 1229, 2967,1*230, 127. * 

1231.315.. 1232.135. . 1233.137. 

22'3*3776*97 ~* " - ' 

321.35. 2. 3 22. 23.0. 323, L3. 7, 324, 31. 6, 32 5, 39.4 

3’2"6»6*27"3",'327 , 6*2. '9 ,*32 8,7 8. "3*i 32*9 , 7 7 .*1, 32 10, 94. 2'"’ "* 

3211. 91. 4. 32 12. 110. . 3213. 106.. 3214. 126.. 3215. 120. 

52*1 6 , 14*27 7 32X7 *, 1 3 4 . , 3'2* L 8' ,"1 5*8”. 73*2X9 , 1'4'9 . ,* 3'2 2*0 , 174 . 

3221. 163. . 3222.190. .3223.177. .3224.206.. 3225.191. 
y27672X2773T2*772'067, 3*22*87278. , 3T297220 ., 3zTo,T54. 

3231.234. . 3232.270. .3233.125. 

41*1, 1.2*1, 4 f2*, 87*29", 413, 9. 2f, 414, 8. 09, 41*5*, *1*8. 6 

416. 21. 4. 417. 21. 4. 418. 22. 0. 419. 22. 0.41 10. 26. 9 

411 1, 26.9, 4112 ,3"1. 0, 4*1X3,31. 0 , 4L*X4, 34. 0, 4~fi5, 3*4.0 

4110.38.9.4117.38.9.4118.42.0. 41 19.42. 0.4120.46.9 
412~1, 46. 9, 4122, 5 1.0, 4123, 5 1.0, 4124, 55. 0,412*5, 55.0 
4126, 58.8,4127,58.8,4128,62.9, 4129,62.9,4130,66.9 
4131*766. 9, 4132.X36. 

421.20.5.422.12.3.423.7.18.424.16.4.425.20.0 
426, 16. 8, 4 2*7, 30787 428, *20. 0,42*9, 3*6. 574210, 23. 2 

42 11, 43. 7,4212, 28. 8, 4213, 52. 7, 4214, 34.9, 4215, 62. 3 

42 1 6. 42. 0. 42 17*, 7T. 6, 4*2X8, 49. 9, 42X9, 8*3. 6, 4220, 577*0 

4221. 96. 6 .4 222. 69. 4. 42 23 . 110.. 4224.8 1 . 6. 4 225, 126. 

42 Si 6 t 96*0 f 42 2 7 1 1 43 • »4220»X15 #j4229» 162 • * 42^0* 136 • 

4231, 183*, 4232, 162., 4233, 103 _ 

5X371773 , 5X2 /l'6 . 4 , 5X3 , 1*8 . 775X4,1 6". 8 , 5 1 5* , *2 6 . 1 
516,22.5,517,32.7, 518,32.7,519,37.2,51 10,39.5 

5 ixrrvf 7T7 51X77737X75X1 3,46. 575T1 47‘fBTl , 5 1 15,49.4 

51 16, 50.9, 51 1_7, 53. l,51_18^51_.5 L 5119,53.7i5120j52_.6 

5X2*1 , 55.4, 5122', 5*3 .*27*5*123, 55 . 6, 5124, 5*2. 9, *5X25 ,54.8 

51 2 6,52.1,5127,53.6,5128,50.6,5129,51.6 , 513 0 , 49.7 

51*31, 48. 7,5132,72.7 
521,10.2,522 ,6.40,523,3.78 
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629,8.83,525,-11.9,526,6.03,527,19.2, 528, 9. 11, 529,29.9 
3210, 11.6, 3211,31.2,5212 ,1 9. 9, 52 13, 38. 3, 3219, 17. 5 
5*2 1"!T, 96. 9, '52*16,21" AT, 3217 ,55.3 15*2*1 8 Afb. 1 , 32*19, 6b • 1 

5220, 2 9. 6, 52 21, 77. 9, 5222, 39. 7, 5223,91.9,5229,91.2 

5225 , 1 OST, 52^6 r9T."6 ;V22 7 ,1 28 . ,'b 2”28 ,'5 0'70~, 5 2 2 9 ,‘l 53'. 

5230.69.3.5231 . 18 3 .. 5232 .82. 9. 5233. 1 08 . 
611,9.51,612,30.9,613,39.9 

619. 29. 2. 615. 9 1.2. 616. 32 .8. 617. 93. 0. 61ft, 39. 2 

6110.90.3.6111.39.5.6112.99.0. 61 13.93.5.6119.96.9 

6115.97.9.6116.99.1.6117.50.5.6118.52.6.6119.51.0 

6120.53.5.6121.52.1.6122 .59.7.6123.52.5. 6129.59.6 

6125. 52.0. 6126.53.8.6127.50.9.6128.52.3. 6129.99.1 

6130, 99. 8, 61 31, 97. 8, 6 13 2 , 96. 7, 51 33, 6ft. 6 

621,7.28,622,9.82,623,2.86 

711.9.67.712.32.0 

721, 7.25,722, 9. 73, 723,2.81 

811.9.67.812.32.0 

1339 , 1.63 

901,3.62,902,9.59,903,5.96,909,5.39,905,7.31 sARcAS FOR 

906.8.29.907.9.16.908.10.1.909.11.0. 910. 11.9 bCONVECTH/E 

911.12.9.912.13.8.913.19.7.919.15.6 J.C0M0UC1 OPS 
921,6. 1 6,9 22,9.02,923,2. 39, 92_9,j> .53 , 9 2 5 , 6. 89. 926 ,9.05 3 A S 
927,1 l7o, 9*2*8 , 5 .09*7*929, 13 A 7930, 6. 12, 931,16.0, 932, 7. 16 TR 0 

933.18. 5. 939. 8.20. 935. 21. 0. 936. 9. 23. 937. 23. 5. 938. 10. 3 6E F U 

939.26.0. 990.11.3.991.28.5.992.12.3.993.31.0.999.13.9 SA 0 K ' 
995,33.5,996,19.9, 997,36.0,996,15.5 ,999 j 3 8.5.950,16.5 IS R C 

951.91.0. 952.17.5.953.7.26 ' 1- E 

19,-1,20,0,21 ,19. , 22,68. ,23,80., 29,22 , 25 , 19,26, 18,27,-6 1 

28,3.2,29, 55, 30*i -8,3 1,5 6, 32, -6, 3 3, 57 ,3 9, -7 

35,-50,36,17 


37,7-66 SNO. UF POINTS FOR TRANSIENT PLOTS 

38,197, 39 , 19,90 , 20,91,21,92, 2 9 ______ 

93,2,99,1, 95, 2'SFOR SKPTST CALLS 

96.0. 97. 1 tFOK STOARY CALLS 

50.0. .51,CUND,52,901 ~ ” 

ENO 


8CO 3ARRAY DATA 

-1, SPACE, 100, ENO 


-2,0.0,1.60,1.90,2.20,2.61 ,2 . 93 , 3. 39 , 3 . 68 ,9.05, 9.93, 9. 8 1 5T0P 
5. 18, 5. 56, 5. 93, 6. 31, 6. 68, 7. 06, 7. 93, 7. 80, ft. 18, 8. 56, 8. 93 »AND 
9.31,9. 68 , 10. 06, 10. 93, 10. 81, 11. 18, 'll". 5 6, 1TT93, 12.31 tBOT- 
12. 68, 12. 92, ENO STOM ROWS OF NODES 


-3, SPACE, 33, END 
-9, SPACE, 33, ENO 

-5, 0.0, 1.60,1. 90, 2‘.20, 2. 61 ,2.82, 2.93,3.09,3.39,3.57 *RDW OF 
3.68,3.79,9.05,9.32,9.93,9.59,9.81,5.07,5.18,5.29, 5.5biNOOES 
5.82,S.93,6.09,6.31~,6.5 7,6.6B,6.79,7.0 6,7.32,7.93,7.59*CON- 
7. BO, 8. 07, 8. 10, 8. 29, 8. 56, 8. 82, B. 93, 9. 09, 9. 31, 9. 57, 9. 684T A I N- 
9.' 79,i0.'o6riO.32,re'.9'3, 10.59,10.81,11.07711.18,11.29 SING 
11.56,11.82,11.93, 12.09, 12.31,12.57, 12.68, 12.79,12.97 SCOQL- 
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-U, SPACE, 100, ENO 

-12. SPACE, 100, END 

-13 

8C0 8RADIAL POSI T10N > I NCHES F-FRDNT ROM, C-OEPTH OF 
■bCDTCDOLAMT T uETEyT^'SrfTRW 

END 

vn 

BCD BTEMPERATUKt DL6RFES F 

TOni 

END 

_ T? ___ - 

8C0 8STEADY TEMPERATURE. I NTE RNALLY-DUOLEO GREGORIAN 
"ffCO - 1.U0E - 0TORT7T FURTH FTDH ' " 

END 

=rf5 

ECO 8AXIAL DISTANCE FROM FRONT SURFACE, I NCHES X-COL. 

BCU 43, Y^CT3C7~ 5, 2~C0XT“2'0 
END 

~Ff~ 3j 

HCD 8T I Mb .HOURS A-NODE 3,K-N0DE 5.C-N0DE 20.D-NODE 1 

TTCO AS sTfe-NCttie" ITS", H-nODE - 1 69'“ ” ~ 

END 

-r»~TO5TTT¥57T37V," _ 3'0GTT5T3'3',‘r: 4~£T6 , 2737) , 371 57E>ID 

-19, .065, . 195,. 3 25 ,.500,. 845,1. A 52, END 

2'07TO"6'5‘,‘: 1957732 5773W. GAO ", 7 3 3 27 1 '.”2 1'2'7L : N D 

-21 

7s“CD — 8T R"A7TSTE'SfT — 'Tt' ^P'ERATURt'TI (iff t’ft N A LUV- COOL E C“G R E GO RI A 
OCO 4N.ULE QUART 2, TURB-.FLOW 

EfTC 

ENU _ _ 

"BCD 3 EXfc'C UT’I ON 

MLTPLY(K2,K3,K10> ^COMPUTE RHO-C PRODUCT 

fiCDAKY ( AT, K1301, K1 302, <1303, Kl'33'4, K1305, <1306, <1307, K1308 
K1309,K1310,K13ll,K13l2,K1313,.<1314,X1315,K1316 

K1 3 17 ,'KT 31 8T<T"3T"9 7Kr3T0TKl 32 l',‘< 1'322', <13 2'3Vra'3'2 4' 

K1325,K1326,<1327,K1328,K1329,<1330,K1331,K1332 _ 

■“ " KL333 1 ----- - 

ARYMPYiKll ,A1,K10,AI) 

uXT<'A‘RYrAl",Ll',:2 _ ,C3,C4,:5;C6,C7,:8,C9;C13,CLL, C12,r.l3,C'l4,Cl'5 
C16, Cl 7 , Cl 8 ,3 19 , C20 , C21 , C22, C23, C24 , C25, C26, C27 , C28 

C 2'9V C 3"0TC 3 1“, C 3‘2 VC3‘3~J ' “ 

BKKAR Y(A1,C34,C35,C36,C37,C38,C3 9,C40,C41,C42,C43._C44, C45 

. ''C4 6';C4'7,'C'48,"C497C50,C5T,'C52,'C53',C54,'C55,C567C5r 

C58,C59,C60,C61,C62,C53,C64,C65,C66) 

BRKAR'Y'CA 1', C 6 7 , C 6 8T C 6 9 ,'C 7 0 ,'C 7 1 , C 7'2, C 7 3 , C 7 4 C 7 5", C 7 6 ,'C 7T , C 7'8 . 

C79,Cao,CBl,C82,C83,C34,CB5,CflS,C87,C68,C89, C90 

C9XrC9'2 _ ,r93;C947C95, , C9'6',‘C977C;3yCW) “ 

MLTPLY(K1334,M0,C331 

OQTAkYCAII K140O , XI 401 ,X1'402 , K1 403 , K1 40 4,< 140'5 , KT40V,K'lVd'7 

KL40B,K1409,K1410,K14U,K141 2,< 1413,XI414,K141S 

" Rl 4 i 67 K1 4 i 7 fK 1' 4 1'87K1 41 9' , XT42"0 , X f42'l , X f42 2 , X'l 4 23“ 
K1424,M42S,X1426,K142 7,K142 8,<1429,K1430,K1431 

KT4T27trr5T37i7r4l4 , k fTSTTfC 1 4 3~£7< 1 4 37 , x 1 4 3 h , x TATV 

K1440, K1441, K 1442, K1 443, K1444,, <1445, X1446.K1 447 
R7l 448", Rl 44 9', X 1 4 5'0 , K 1 45'1 , K 1 45 2 , X 1 453% K 1 45 4', K 1 455 
M456,M457,X145B,M4S9,K1460, X146 1,X1462,K1463 
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K1464,Kl465,K1466,K1457,K146 8,ia469,K1470,K.14/l 
K14 72,Kt473,K1474,KL475,K1476,lU477,!U47a,K1479 
K1480.K1 481,l<140'z , iaT3T,K14f4 1 4615 14B6,K 1487' 

K14 88 , K1489 ,111490) 

ARYMPY(Rl2,Ai»"kT0,AT)' 

BRKARY( A1 , C100 » C 101 , C 102 , 003 , Cl 04, C 105, C106, C107, CIOS 

Cl o 9 , c 1 1 o , c i fi , c 1 T 2 , c l iT,t 1 1 4 , c ITS', C i 1 6 ',t 1 1 ? 

Cl 18, Cl 19, 0 2 0, 02 1,0 22, 0 2 3, 024,0 25,026 

C12TTC128 ,029 ,Clj6,(7l'3l , (Tl 5 - 2 , £T 133 , - Ci 34, Cl 35 
C136,C137,C138,C139,C140,C141,C142,C143,C144 
04$, 040 047“, 048, 0 49, 05 0,051, 052, 053 
C154, C 1 58 1 Cl 56 , C 15 7 , Cl 58, C159, C 160 , C161, C162 

C 16 3 , C 1 r4 , C 1 6'S7C1 6 6 , 1 11.7 , ere 8 , C 1 69 , C 1 7 0 , Ei 7 i 

C172,C173,C174,C175,C176,C177,C178,C179,C180 
Cl bT, £182,08 3, CT84, Cl 85, C1B 6, Cl 87, 088, 08$ 

Cl 90) 

bCDAR Y ("A I , K1 1 r, K 1 f2 , K 1 1 3 , K 1 1"4 ,'lCl 15 , _ K 1 1 6 , X.Tf7 , K 1 1 8 , K 1 1 9 

K1110,Kllll ,M 112, Kill 3,K1114 t K l 11 5 ,_K1 1 1 6 

KH17,K1118,K1119,KH20,K1121,<ri2^,KlT5T 
K1124,K1125,K1126,K1127,Kll28»i<1129,K1130 
K1131 , K1 132 ) 

ARYMPV(K13,A1,K1,A1I 

BRKARYIAl.GUi ,G112,G113,GU4,G1 15, Gil 6, G1 17TG118, G1 19 
Gill 0, G1 111,3 1112, Gill 3, SI 114, 31115,61116 
GTl lT7G 1118,3 1 119,01 120, Gfl21,Gfi2'2,Gl 123 
G1124,G1 125,G1 126,G1 127 ,G1 12 8 , G1 129, G1 130 
G1131 , G1 132 ) 

BLDARY(Al',K121,K122iK123,K124,K125,)U26,A127,K128,Kl29 
K1210,K1211,K1212,K1213,K1214, A1215.K1216 

K1217,K12I8,K1219,K1220,K1 221, <1222, K1223 

K1224,K1225»K1226,K1227,R1228, <1229,0 230 

K1231.K1232.Q233) 

ARYMPYt Kll i A1 ,K1 , A 1 ) 

8RKAKY(Al,G121,G122iG123,G124,G125,G126,G127,G128,G129 

G1210,G1211,G1212,G1213,G1214,G1215,G1216 

G1217,G1218,G1219,G1220,G1221,G1222,G1223 

G12Z4,G1226,G1226, G1 22 7 , G1 228 , 31229 , G1230 

G1231.G1232.S1233) 

MLTPLY(K2233,K1, 32233) 

8L.0A8VI A1 ,K321,K322,K323,K324,K3 25,K326, <327 , K32B, K329 
K3 21 0, K3 21 1,K3212,K3 213, K3 21 4,. <3215,63216 

K32I7,K3218,K3219,K3220,K3221,K3222,K3223 

R3224,K3225,K3226,K322 7 , K.3 22 8 , <3229, <3230 
K3231.K3232.K323J) 

ARYHPY(K11,A1,K1,A1) 

BRKARYtAl, G3 21, G322, G323, G32 4, G3 25, G 326,0327, G32B ■ G329 
G3210,G3211,33212,G3213,G3214,G3215,G3216 

G3217,G3218,S3219,G3220,G3221,G3222,G3223 

G3224,G3225,G3226,G3227,G3228,G3229,G3230 

G3 231, G3 232,33233) . 

MLTPLY<K621,K1,G621) 

HLTPLY(K622,K1,G622) 

HLTPLY ( K623 , K 1 ,G623) 

HLTPLY(K711,K1,G711) 

MLTPLY(K712,(U,G712> 

MLTPLYIK721, 0,6721) 
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MLTPlYIK722,K1,G722> 

HLTPLY<K723,K1,G7231 

MLTPLY(K8U.K1,GB11) 

MLTPLY<KB12,K1,G812> 

BLDARY ( Al, <411, <41 2, K41 3, <414,4415, <416, 4417, <41 8, <419 

<4110, <4111, <4112, <4113, <4114, 44 115, 44116 

<41 1 7', K4TfS,<4U9, <412 0, <412 1,44122, <4123 

<4 124, <4125, <4126, <4 12 7, <4 1 2 8, <4129, <4 130 

K4131,K4132> 

ARYMPY(K13,A1,<1 ,A1 ) 

BRKARYf A1,G411,G412,G413 ,G41 4, 34 15, G416 , 3417, G4 1 8.G419 
G4110,G4in,G4112,G4113,G4114,G4115,G4116 
G'4Tl 7, GUO), 54 119, G4120.G41Z1, 54122, G4123 

G4124, G412S, 34126, G4127, 04128,54 129, G4 130 

G4131.G4 L32) 

BLOARY ( A1 , K421 ,<422, K423 , <424, <425, <426,. <427, <428 ,<429 
iv42T0 , <42 fl, <4212, <4 213, <4214, <4215, <4216 
<42 1 7, <4 21 B, <4219, <4220, <4221, <4222* <4223 

<4224, <4225, <4 226, <422 7, <4228,44229, <4230 

K 4231, <4232, K4233) 

ARYMPY(~KTT7A1 , <1 , A1 ) 

BRKARY(Al,G421,G422> G423 ,G424, G4 25, G426,S427 , G428 ,G4?9 
|j4"2"1 0 1 G42TT7G 4212, G 4 21 3', 0^2147 54 2 15 , G 4 2 1 6 
G4217,G421B, 54219, G4220,G4221, 34222, G4223 
G4224.G4 225734226, G422 7, G4228.G4229, G4230 
G4231.G4232.G42331 

— BTC5SRYTAT7K5 1 1 , < SlT, <513, <51 4, <5*15, <51*6 , 45 1 7 , K 5 IS , KSl'9 
<5 110, K5111, <5112, K5113, <5114,45115, <5 116 

< 5 Ti 7 , < rrrerirs 1 19, <5120, < 5 1 2 1 , < 5 1 2 2 , <s 12 3 

<5124, <5125, <5 126, <5127, <5128,45 129, <5 130 

K5T3TTR rrsTi 

ARYMPY(<13»A1,<1,A1> 

^7fKXTYTST7(J5Tr7G'5T2T5'5r3TG'S'r4TG5T5TG"5'r<j'r5'5'17“, G5TB ,‘G5T9 
G5UO,G5111 ,G5 112,55 113, G5114, 55 tl5,G5 116 

35 rrr 7G‘5Tra',Ts‘iT9 ; g'syzs rt n? rv 5 512 2 , g 51 2 3 • 

G5124.G5 125, 55 126, G5127, 55128,35 129, G5 130 
GvmTcvn^i 

BLnARY(Al,<52t,<S22,K523,<524,K525,<526,<527,K528,<529_ 
<5216", K5'2'l"ir<52'f2',K52'13',K5214, <5215, <5216 
<52 17, <5218, <521 9, <5220, <522 1,45222, <5223 
R5'2T4,<5'2'2'5,<5'22'6,<5227',<5228, <'5229, <5230 
<5231, <5232, <52331 

"SRYRPTrKTI 7 AX,ia 7 AT) 

BR<ARY(A1,G521,G522,G523,G524,G525,G526,3527,G528,G529 

G52Y0TG5 2T1T552T.2, G521 3 ,G5”21'4,G52 15i G5216 ’ ’ 

G5217, 65 21B, 35219, G5220,G5221, 55222, G5223 

C 5 2Y47GT22 575 5 2 2~6 G 522' 7,'C52 2 a , 3 5 2 2 9 , G 5 2 3 0 ' 

G5231, 05232,35233) 

_ ii'CTOTvr)riT<6'nTK:^'2Tr6T3TiC6Yvr<'5T5;<6r6 _ ,X6'i7', <6ie 

<6 1 10, <6111,461 12, <6 11 3, <611 4,46115, <6 116 
<Sl 1 7, Kisri 8, 461 19; <6120, <61.2!, 46122 ,<6123 
<6 124, <612 5, 46126, <6 127, <6 128, 46129, <6130 
<6131 ,<6132 } 

ARYMPY(K15,A1,K1,A1) 

"TTR'KTrR.'Y'nri'TGSYlVGS T2TG61T, Cbl 47 G6'15 ,'G61 6 , 3617 , G6 1'8 

C6110, 06111, 36112, Gbll 3, G6114, 36115,06116 
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G6117,G6118, 36119, Gbl 20 , G612 1, G6122, G6123 
G6124.G6125,G6126.G6127,G6128,G61 29,G6130 
G6"131,G£>132) ‘ ' 

BLDARYI A1,K901,K902,K903,K904,K9D5,K906,X907,K908 
^TTOTfCsTO , K9lT ,'R9'l 2 , K9 13 , K9~T4) 

ARYMPY(K16,A1,K6,A1> _ 

ARYO I VTI< 1 6,Al,Kt , A1 > 

BRKARYI A1.G901.G902,G903,G904,G90S.G906,S907,G908 

C9WTG9'1 0, (J9T17GT1 2 , G9 13'V G9 1 4 ) 

CSGDMP 



QPNPL7 

CTflD'S'S 

BLDARYI A3, 71, 72, 73, 74, 75, 76, 77, 78, 79, 710, 7 11, 712, 713 _ 

— TlATrfSTT r67rT77T f & , r 1 9 ,1 2 o , t 2 1 , t i z , r 2 3 , r 2 4 

T25,726,727,728,T29,730,731,732,733) 

BLDARY(A4;'T183,tl84,TlB5,T152,rL53 r ri54, ri56,T157,T158 

T159, 7160, 7161, 7162,71 63, T164, 7165,7166,1 167 

rf5"8T7169,717 0,T171,7172»Tl7 3,T174,7175,7176 

71 7 7 , T 178. T 179,7180.1181.1182) 

ET'DATTyI A 6 ,T 1“ 00 ,7101, n.02,7103, II D4, 7194,. 72 4 7 , 7 192 ,'7 105' 

T198, 7248, 719b, 710b, 72 02, 7249, 7200, 7107, 7206 

7250,7204, 7108,7210, 7251, 7208,7109, 7214,1 252 

T2 1Z, 71 1 0,721 8,7253, 72 16, T 11 1,T222, 7254, 7220 

T112, 7226, 7255, 7224, 71 13,723 0,1 256, 7228, 7114 
T234, 7257, 7232, 7115, 7238,7258,7236,7116,7242 
T259T7240, 711 7, 7246, 7260,7244,71181 
PLOTX2(K19,K20,K21,K22,K23,K24,A13,A14,A15,K11,A2,A3) 
RL0TX2(K20,K20,K21,K22,KZ3,K26, A13.A14 ,A15,KL1, A2, A4 ) 
PL07X2 (K20>K20,K21 ,K22 , K 2 3 , K25 , A13 , A14 , A15,K2 7 , A5,A6) 
BLDARV'i"AT,T3T736,769,7102,T121,7185,ri88,7155T 

BLDARYI A4, 75,738, 771, 7104, 7123, 7 153) 

BLDTR'Y’rA'6T721T,T5ljT7’8'6T7 2 1 9 , T22 1 , T ll 8 , 7 1 6 9 ) 

PL07X2 (K19,K20,K2B ,K22, K23 ,K29,A16,A14,A15,K30,A18,A3) 

pL07X2IK20,K26,lC2'8";K22,R2i,R3l,4L6,Al4,Al5,K32,Al9,A4r 

PL07X2(K20,K20,K28;K22,K23,K33, A16,A14, A15,K 34, A20,A6) 

SYf s q £“ I K4 , kTTTt 1 > 

S7FSEP I K197,7IHBM0) 

S7FSEPIK19a,OU7PU7) 

S7FSEP IK199,CSGFAC ) 

CNFRWD 

PL07X2IK19,K20,K35,K22,K23,X36,A17,A14,A21,K37,A1,A7) 
PL0TX2(K26,ir2O",'K35',K22,K23,X38,A17,A14,A21,K37,Al, AB ) 
PL07X2(K20rK20,K35,K22,K23,K39,A17,A14, A21,K37,A1, A9 ) 
PL07X2(K20,K20,K35,K22,K23,K40, A17, A14 , A21, K37, Al, AllO 
PL0TX2(K20,K20,K35,K22,K23,K41, A17.A14, A21,K37,A1, All) 
Pl67X2(K20,K20,K35,K22,K23,K42,A17,A14,A21,K37,A1,A12) 

Pl7NP 

SKP7S7IK20,K45,K44,K44) 

ENDFIL 

END 

BCD 3VAR1ABLES 1 

BLDARYI A1,K925,K926,R927,K928,X9?9 

K930.K931.R932 , K933,K9 34,X93 5 , K936, K937 , K938 

X939,K940,K941,K942,K943,X944,X945,K94b,K947 

X948,K949,K950,K951,K952,K953) 



Table 13. Chrysler Improved Numerical Differencing Analyzer (Contd) 


ARYMPYOUB.Al ,K5,A1) 

BRKARY( Al, 05,06, 07, 08, 09, 010,011 ,012.013, 014,015,016, Q 17 
018,019,020,021,022,0231024.025. 026,027,028, Q29 

030,031, 032,033) 

ADD 105,06,07,08, Q9.Q10, OH, 012, 013, 014,015,016, Q17,Q18 

01 9, 020,021, 022, 023, Q24.Q2S, 026, 027, 028, 029,030, 031 

032,033, *50) 

BNO 

BCD 3VARI ABLES 2 

END 

BCD 3CIUTPUT CALLS 

INDEX1K46 , K47 ) 

ST0ARY(K46,T3,A7) 

ST0ARY1K46,TS,AB) 

STOAR Y ( K4A ,T20,A9) 

ST0ARY(K46,T155,A10) 

ST0ARY(K46,T153,A11) 

STOARYt K46 ,T169,A12I 

STOARY IK46 , TI MEN , A1 ) 

i-ND f£3~7 ) 
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program the correct value will be calculated and loaded into these locations, 
except for those nodes that are either a boundary node, which is noted by a 
minus sign in front of the node number, or an arithmetic node which is noted by 
a negative capacitance value. 

The next block of data is computer generated and assigns numbers to each 
node listing first the diffusion nodes (those nodes with capacitance) then the 
arithmetic nodes and last the boundary nodes. 

Following this block of data is the "BCD 3 CONDUCTOR DATA". In this 
block the first value is the conductor number, followed by paired numbers which 
are node numbers which have this conductor value between them. The final 
number lists the value of the conductor. Again the conductor value input is zero. 
The correct value will be calculated at a later time and will replace the zero. 

The next block of data is the "BCD 3CONSTANTS DATA". This block of 
data lists the constants which are necessary to run the program, geometrical 
factors, and physical characteristics of the materials being considered. It is 
from this stored constant data that the thermal capacitance and thermal con- 
ductivity are calculated. 

Following the "BCD 3CONSTANTS DATA" block is the "BCD 3 ARRAY 
DATA", Arrays are established which allow operations on groups of data, allow 
storage area for data to be plotted and store labels for tables and plots. 

The next block of data is the "BCD 3EXECUTION". It is in this block that 
the many subroutines are calLed up to do the computations to completely define 
the network and finally call upon the thermal analysis subroutine which is 
desired. In the sample case the card CINDSS calls upon a steady state sub-' 
routine to get the steady state solution. Following this is a card GNFRWD 
which calls up a transient subroutine. 

Other instructions are given in this block also. The computer is 

/ 

instructed to print out certain data values desired where those values are con- 
stant. The computer is also instructed regarding plots for the steady state, 
and the transient solutions. 

Following the "BCD 3EXECUTION" block is the "BCD 3 VARIABLES 1" 
where it is required that the data of impressed source (solar, electrical, etc. ) 
for any node be located. The impressed source is zeroed out after each iter- 
ation and reloaded at the beginning of each iteration. The impressed source 
data can either be constant values, time variant or temperature variant. If 
a print statement is included in this block, the command is honored after each 
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iteration and a printout will appear. This can be utilized to good advantage 
when a problem in the program exists. Since the printout appears after each 
iteration, much paper will be used so it is not recommended that this be done 
except for trouble shooting purposes. 

In the "BCD 3 VARIABLES 2" block post solution operations are per- 
formed. Examples of such operations are the rate of heat flow between nodes 
and the integral of this heat flow rate over a given period of time. Any calcu- 
lations which are to be done after an iteration has been completed is to be 
called out in this block. 

The final block in the program is the "BCD 30UTPUT CALLS". Through- 
out the program much data is stored in many locations, but it is not printed out 
or plotted unless a specific instruction is given. In this particular program 
data are stored by the command "STOARY" which the plot subroutine in the 
"EXECUTION" bLock utilizes for the transient plots. 

It was pointed out in the "VARIABLES 1" bLock that a print statement 
there would cause a printout after each iteration. A print statement in the 
"OUTPUT CALLS" block causes the data to be printed out after each time step. 
Each block must end with an "END" card. 

The output of the program can be either in printed or graphical form. 

Table 14 shows the temperature listing for each node in the printed form. The 
first three columns indicate the relative node numbers of the temperatures 
listed in that line. For instance, the first line reads 1 THRU 5, The five 
numbers following are the temperatures of nodes 1, 2, 3, 4 and 5. 

■ Figure 18 is a hand plot of the temperature of the front and back surface 
of the plane mirror. The back surface curve is a constant 70 °F. There are 
three curves for three different materials - Quartz, Invar and ALuminum. This 
is the thermal data which was given to the stress analysis group for inputs to 
their program. Figure 19 is an example of a computer -generated plot. 

Considerations other than thermal governed the configuration of the mir- 
ror, In the early stages, two configurations were considered. The plane mir - 
ror, Figure 17, which has been discussed before, and a tapered mirror shown 
in Figure 20. Figures 17 and 20 show a cross-section of one half of the mirror. 
The rectangles shown in these figures are revolved about the centerLine to gen- 
erate a quadrilateral torus. Figure 21 is a top view of the mirror showing the 
surface of revolution and node boundaries. 
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Table 14. Overall Analysis 
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Figure 20. Tapered Gregorian Primary Mirror 






















DOTTED LINES REPRESENT INTERNAL NODE BOUNDARIES 
CROSS-HATCHED AREA IS SHADED BY SECONDARY 
RADII TO OUTER BOUNDARIES OF NODES: 


r 70 = 1/2 in. 

ig = 4.71 in. 

Xy = 11.0 in. 

= 14.9 in 

r 8(| = 1.5 in. 

— 6.28 in. 

Tg = 12.6 in. 

i-j ^ = 15.7 in 

=2.36 in. 

= 7.85 in. 

r 9 = 14.1 in. 


r« • = 3. 14 in. 

= 9.42 in. 




Figure 21. Radial Subdivision of Plane Gregorian Mirror 
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The diameter of the tapered mirror was the same as the plane (80 cm) 
and tapered from 1-1/2 inches thick at the center to 1/2 inch at the periphery. 

The boundary conditions for the thermal analysis are as follows: 

1. One solar constant was uniformly incident on the front surface with 
the exception of an on-center spot 12 cm in diameter. This spot 
represents the effects of shading by the secondary assembly, 
excepting the web. 

2. A solar absorptance of 0, 14 was used for all cases. 

3. Heat was removed only at the bottom surface of the mirror. (The 
method of accomplishing this boundary condition is not defined). 

Four different temperature distributions were specified for the bottom 
surface: 

1. A uniform temperature along the entire radius, equal in value to 
the initial (uniform) temperature of the mirror. 

2. A uniform temperature along the radius from the edge of the shaded 
spot to the periphery, equal in value to the initial mirror 
temperature. 

3. Temperature increasing parabolically from the initial mirror tem- 
perature at the periphery to 3°F above this temperature at the edge 
of the shaded spot. 

4. Temperature increasing parabolically from the initial mirror tem- 
perature at the edge of the shaded spot to 3°F above this tempera- 
ture at the periphery. 

The results of this investigation are shown on plots in Figures 22 to 39. The 
four temperature distributions listed above are illustrated by the lower curves 
in Figures 22, 24, 26, and 28 respectively. 

It should be noted that a value of 70 °F was used as the initial value for the 
plane mirror, while 50°F was used for the tapered mirror. This difference 
has no significance; and it does not affect the thermal distortions, since they 
depend on temperature differences rather than the scale employed. 

The results for each combination of geometry and back surface tempera- 
ture distribution are presented by a pair of figures. The first figure of each 
pair presents the steady radial temperature distribution in terms of the tem- 
peratures of the front and back rows of nodes. The top curve pertains to the, 
front row and the bottom curve to the back row, in every case. No attempt 
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Figure 25. 
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Figure 32. 
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was made to obtain smooth curves for the steady state results, so 
unrealisticaLLy sharp corners appear in these curves. The second figure in 
each pair presents the transient results in terms of the temperatures of three 
nodes in the top row. The curve drawn through the +'s is for the center node; 
that through the 0's is for the node just inside the edge of the shaded spot in the 
case of the plane mirror, and for the node centered about 2 inches beyond the 
edge of the shaded spot in the case of the tapered mirror; that through the *'s 
is for the node at the periphery of the mirror (see Figure 23 for example). 

It is noted that the temperature difference between the front and back 
rows of nodes of the plane mirror approaches .a nearly constant value of about 
5, 3°F for radii greater than 4 inches, irrespective of the temperature distribu- 
tion along the back surface. This indicates that radial conduction is not a very 
important factor in establishing the temperature distribution, and this indica- 
tion is supported by other features of the results for both the plane and tapered 
mirrors. One feature which supports it is that the steady temperature of the 
shaded region is much lower if this region is cooled from the back than if it is 
not cooled. 

The transient plots for the plane mirror indicate that the temperature 
near the periphery approaches steady state closely at the end of 0. 3 hour, 
irrespective of the temperature distribution on the back surface. The temper- 
ature near the edge of the shaded spot approaches equilibrium less rapidly; 
because it depends on radial conduction; and it is more sensitive to the temper- 
ature distribution on the back surface. The temperature at the center of the 
mirror approaches equilibrium very slowly because of the low thermal diffus- 
ivity of the material. 

The results for the tapered mirror (Figures 30 to 39) follow the same 
trends as those for the plane mirror, although the distributions are somewhat 
more complicated because of the radial variation in mirror thickness. The 
radial distribution of the temperature difference between the front and back 
rows of nodes is virtually independent of the radial distribution of back surface 
temperature. The plots of transient temperatures show that equilibrium at the 
three radial locations is approached in the same order as for the plane mirror. 
However, the node near the periphery of the tapered mirror approaches equi- 
librium more rapidly than the corresponding node of the plane mirror, because 
of the reduced thickness at that position, while the node near the edge of the 
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shaded spot approaches equilibrium less rapidly, because of the increased 
thickness at that position. 

After the 80 cm primary mirror investigation was completed, the mirror 

size was reduced to 65 cm diameter. It was recognized that the ATM con- 
BXTJ 

straint of 0. 6 ■ heat transfer to the spar for each attachment point would 

require that the mirror temperature would have to be held close to the 70 °F 
ATM spar temperature. Also, the problem of getting a system to accomplish 
the controlled temperature on the back surface seemed insurmountable. 

For these reasons an investigation of a 65 cm diameter mirror with inter- 
nal cooling channels for the purpose of circulating a cooling fluid was begun. 

The material U. L. E. (ultra low expansion) fused silica made such a design 
possible due to its fusibility. Figure 40 shows a cross section of one half of 
the mirror. 

A detailed analysis was made of the steady and transient temperature dis- 
tributions in an annular segment of the 65 cm Gregorian primary with an inter- 
nal cooling channel. The purpose of this analysis was to obtain a detailed pic- 
ture of the temperature distribution around the cooling channel, so that the effect 
of local thermal distortions on the mirror figure might be better assessed. 

The configuration of a cross section of one -half the mirror, shown in 
Figure 41, is the same in ail essential thermal respects to the latest version, 
namely, it is 65 cm in diam. and tapers from a thickness of 2 in. at the center 
to 1 in. at the periphery, A circular cooling channel of l/4-in. diam. spirals 
from the edge of the shaded spot to the outer edge, with successive turns 
spaced 3/4 in. apart on centers. The tube centers are 1/2 in. from the illu- 
minated surface of the mirror. 

The segment of the mirror selected for detailed analysis was the thickest 
section of the illuminated portion of the mirror, as shown between sections A -A 
and B-B in Figure 41. The selection of this segment was made to obtain 
results at a position where local thermal distortions should be greatest. The 
spiral nature of the cooling channel was ignored so that the segment could be 
treated as circular in shape for simplicity of analysis. The resulting segment 
is in the form of a ring with a concentric circular channel, as shown by the 
sectional view in Figure 42. To further simplify the analysis, large-scale 
radial temperature gradient's were ignored, which yielded symmetry of the 
temperature about the cylindrical surface passing vertically through the center 
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Figure 41. Set-up for Detailed Analysis of Tapered Gregorian Primary 
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Figure 42. Section of 65cm, Tapered Gregorian Primary 
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of the tube, so that only half of the segment need be analyzed. The segment was 
subdivided vertically and horizontally into 76 elements, as shown in Figure 42, 
Convective resistance between the wall of the channel and the coolant fluid was 
taken into account, since a rough analysis indicated that this resistance was 
significant compared to that due to conduction through the mirror. Results 
were obtained for both laminar and turbulent convective resistances, (Laminar 
flow corresponds roughly to a 3°F rise in coolant fluid temperature between 
inlet and outlet to the primary mirror, whereas turbulent flow corresponds to a 
rise no greater than 1. 5°F. ) 

The coolant fluid is represented by node 99 in Figure 42. Although the 
temperature of the coolant increases as it flows through the channel, the varia- 
tion in one turn is small, so it was ignored. To help account for the transient 
thermal process undergone by the entire primary mirror cooling system at the 
start of viewing, the temperature of the coolant was varied with time in some 
cases. The types of variation will be described later in discussing the results. 

Results are given in Figures 43 through 45 for the temperature distribu- 
tion in the mirror segment due to solar irradiation and turbulent flow of a cool- 
ant fluid having a steady 70 °F temperature, corresponding to an infinite flow 
rate. Figure 43 shows results for the steady radial temperature distribution 
between adjacent turns of the cooling channel. The upper curve in Figure 43 
gives the temperature variation of the top row of nodes, from node 1 on the left 
to node 5 on the right, The lower curve in this figure gives the temperature 
variation of the row of nodes directly opposite the tube, from node 37 on the 
left to node 40 on the right. It is observed that the radial variation in tempera- 
ture between adjacent turns of the channel is rather weak, particularly near the 
illuminated surface. This observation applies for all the steady results pre- 
sented herein. 

The vertical variation in the steady temperature distribution is shown in 
Figure 44. The upper curve presents the results for the row of nodes located 
approximately midway between adjacent turns, extending from node 5 at the 
top to node 80 at the bottom. The lower curve presents the results for the row 
of nodes extending from node 1 at the top to node 76 at the bottom. The two 
curves show that the temperature drop from the front surface of the mirror to 
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the back is about 3°F. The variation is smooth except for the portion of the 
bottom curve which accounts for the presence of the channel. An interesting but 
unimportant feature of this curve is that the temperature rises slightly from the 
wall of the channel to the back surface of the mirror. In other words, not all of 
the solar heat absorbed at the front surface of the mirror is transferred to the 
coolant fluid before it is conducted to the region below the tube. 

Figure 45 presents the transient temperature variation at four locations 
in the mirror. The temperature is initially 70°F throughout. It is observed 
that at nodes 1, 5, and 40 the temperatures have very nearly reached steady 
state (indicated by corresponding symbols on vertical axis) about five minutes 
after the solar flux has been suddenly applied. The temperature at node 80 
approaches steady state less rapidly, but this is relatively unimportant because 
of the small change in temperature at this node. 

Figures 46, 47, and 48 present results of the same type as those just dis- 
cussed, except that the convective resistances are typical of laminar, rather 
than turbulent, flow. The higher steady temperatures shown in Figures 46 and 
47, compared to those in Figures 43 and 44, are due to the larger temperature 
differences between the channel wall and the fluid required to transfer the 
absorbed solar energy to the fluid. In other words, for a given temperature 
distribution prior to solar irradiation, the temperature at any point in the solid 
rises to a higher level with laminar flow through the coolant channel than with 
turbulent flow. The transient results for laminar flow, presented in Figure 48, 
indicate that steady state is approached much less rapidly than with turbulent 
flow. This is due to the higher resistance to convective heat transfer which 
occurs with laminar flow, compared to turbulent flow. Not only does the tem- 
perature at any point in the mirror rise higher with laminar flow, but it does so 
more slowly. 

In actual operation, the entire coolant system, including the radiator, 
pump, accumulator, tubing, etc. will undergo a transient thermal process after 
a change from nighttime to daytime conditions. Hence the coolant fluid will not 
remain at a fixed value, as assumed in obtaining the above results, but it will 
go through a transient also. To obtain a reasonable estimate of the variation of 
coolant temperature with time would require thermal analysis of the entire 
system. This is beyond the scope of the present investigation. To gain some 
feeling for the effects of transients in the coolant temperature, two types of 
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Figure 47. 
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transient variation were assumed. One was a step increase of 2°F in the 
coolant temperature at the same instant as the start of solar irradiation. This 
type of variation gives a lower limit on the time required to reach steady state. 
The other type of variation was a rising exponential, given by 

T = T + AT fl - e —t / tc:~j (j) 

qq 0,qq L J 

(when t = t^; - e = .632 = 63. 2%) 

where T is the initial coolant temperature, AT is the maximum rise in 
o, qq 

temperature above the initial value, t is time, and t is the time constant for 
the transient change, i. e. , the time required for 63. 2% of the change to occur. 
The rising exponential is a more realistic type of variation since it predicts an 
asymptotic approach to steady state. By varying t , one may alter the rate at 
which steady state is approached. In the results presented here a value for 
AT of 2°F was used in equation (1), to permit direct comparison with the results 
for the step change in coolant temperature. A value of 0. 1 hr. was used 
throughout for t . 

The steady state results in Figures 49, 50, 51 and 52 are all 2°F higher 
than the corresponding results presented above for a steady coolant tempera- 
ture. This follows logically from the fact that the steady fluid temperature is 
2°F higher than formerly. The transient results for a 2°F step rise in coolant 
temperature with turbulent flow, presented in Figure 53, indicate that steady 
state is nearly achieved at the end of about five minutes, as before. When the 
coolant temperature rises exponentially, the mirror temperatures approach 
steady state less rapidly, as shown by the results in Figure 54. This is pri- 
marily due to the continual change in coolant temperature during the period 
shown. 

The transient results for laminar flow are more sensitive to transient 
changes in coolant temperature, as demonstrated by the curves in Figures 55 
and 56, Both a step change and an exponential change in coolant temperature 
cause a significant increase in the departure of the transient from the steady 
state values at the end of 0. 3 of an hour. 
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Figure 49. 
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Figure 53. 
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The transient and steady state temperature distributions for the entire 
65 cm primary mirror made of U. L. E. with internal cooling channels are 
shown in Figures 57 thru 74. The results of this study are complementary to the 
detailed analysis which describes in detail the temperature distribution in the 
vicinity of one cooling channel. As in the detailed analysis, the convective 
resistance between the channel wall and the bulk of the coolant fluid is taken 
into account in order to gain a more realistic picture of the transient and steady 
state temperature distributions. 

The major assumptions on which this analysis is based are as follows: 

1, The cooling channels are concentric, rather than spiral. This 
corresponds to zero spacing between infinitesimal channels and was 
used to permit the temperature and thermal distortion distributions 
to be treated as axially symmetric. Due to the limitation in the 
number of nodes which can be handled by the thermal distortion 
program, the analysis would have been impossible without this 
assumption. It would have been possible to thermally analyze the 
nonaxially symmetric situation. However, the time required to set 
up the data and carry out the computations would have been much 
longer, because several times as many nodes would have been 
required. 

2. The fluid temperature was assumed to be either uniform in the 
radial direction or to vary in a continuous fashion. Continuous 
variation is consistent with a spiral cooling channel pattern. It 
might be impossible to achieve it with concentric channels, but this 
is thought to be unimportant, since the desire was to simulate the 
spiral pattern as closely as possible. Two directions of flow were 
considered. If the inlet is at the inside, the fluid temperature 
increases parabolically from inside to outside, as shown by curve 1 
in the sketch below. If the inlet is at the outside, the temperature 
again varies parabolically, but as shown by curve 2. Since laminar 
and turbulent flow correspond to different flow rates, the difference 
between inlet and outlet temperatures is not the same for the two 
types of flow. A 1°F difference was used for turbulent flow, which 
difference corresponds to a Reynolds number of 6380. A 3°F dif- 
ference was used for laminar flow, corresponding to a Reynolds 


- 115 - 



750-7 



number of 1600. These Reynolds numbers pertain to a fluid of 
80/20 methanol/water and a circular cross section flow channel, 

1/4 inch in diameter. 

3. The cooling channel was assumed to extend from the periphery of 
the central shaded spot to the periphery of the mirror. This is not 
consistent with the latest version of the primary mirror, which 
shows the spiral extending into the outer portion of the shaded spot. 
However, the thermal distortion program will not permit treatment 
of more turns of the channel than the number considered, namely 14. 
A note of caution is in order, however, since the mirror is sup- ’ 
ported at the hub, and this portion comes to equilibrium quite 
slowly, in a relative sense, as shown by the attached results. In 
other words, cooling the shaded portion might decrease the time 
required for the mirror to approach an equilibrium position, 
although the distorted figure may be insensitive to whether or not 
this portion is cooled. 

4. The effects of an internal radial inlet or return channel in the pri- 
mary mirror were not considered, since this channel would destroy 
axial symmetry. The necessity for assuming such symmetry has 
already been explained. 
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5. The cooLing channels were spaced 3/4 inch apart on centers. To 
simplify the thermal distortion analysis, the 1/4-in. diam. circu- 
lar channels were replaced by square channels of equal cross- 
sectional area. This may have a mild effect on the temperature 
distribution, but any such effect is thought to be relatively 
unimportant. 

6. Heat is removed from the mirror solely by the coolant fluid. That 
is, losses through the mounts and by radiation from the surfaces of 
the mirror were neglected. It is thought that this assumption has 
negligible effect on the temperature distribution for the present 
material and configuration. However, this does not necessarily 
hold if conductive mounts are used to support the mirror at its hub 
and if the material of the mirror is changed to one which has a 
higher thermal diffusivity. 

The subdivision of the mirror for the thermal analysis was essentially 
the same as for the thermal distortion analysis. However, a different node 
numbering scheme was used. Figure 40 shows the discretized mirror and num- 
bering scheme. (It may be noted that the notches in the hub were disregarded 
in the thermal analysis. The effect of this assumption on the results of the 
analysis is thought to be negligible. ) Nodes with zero thermal capacity were 
placed at the four surfaces of each square flow channel to permit more accu- 
rate and more detailed results to be obtained. The temperature of the fluid in 
each channel was represented by a separate node to permit specification of the 
fluid temperature distribution along the radius of the mirror. 

Representative results are presented in Figures 57 through 74. Each group 
of three figures (Figures 57-59, 60-62, 63-65, 66-68, 69-71, 72-74, presents 
the results for one type of flow and one fluid temperature distribution. The 
type of flow specifies the magnitude of the convective resistance, whereas the 
fluid temperature distribution establishes the base for the temperatures in the 
mirror. The first figure in each triplet presents the radial steady-state tem- 
perature distribution. The upper sawtooth curve, denoted by symbol F, gives 
the distribution through the upper, or front, row of nodes from node 1 through 
node 33. The curve denoted by symbol B gives the distribution through the back 
row of nodes from node 183 through node 182. The lower sawtooth curve, 
denoted by symbol C, gives the temperature distribution through the row of 
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nodes in line with the cooling tubes extending from node 100 through node 118. 
The relative maxima in curve F correspond to locations midway between tubes, 
whereas the relative minima correspond to locations centered (radially) above 
the tubes. Thus, the temperature difference between the maxima and minima 
gives approximately the fluctuation in surface temperature due to the presence 
of the internal cooling channels. However, the magnitudes obtained for these 
fluctuations from the present analysis are approximately twice as large in ail 
cases as the fluctuations obtained from the detailed analysis. The results of 
the detailed analysis are probably more accurate due to the finer subdivisions 
used therein. In any event, the maximum fluctuation given by the present analy- 
sis is about 0, 25 °F and occurs for turbulent flow. The relative maxima in the 
curve denoted C correspond to locations midway between tubes, whereas the 
relative minima are fluid temperatures. Thus a curve drawn through the 
minima would represent the fluid temperature distribution with radial position 
on the mirror. Such a curve was not drawn, to avoid confusion. The points 
denoted C which are intermediate to the maxima and minima represent the tem- 
peratures of the coolant channel walls. The temperature differences between 
the channel walls and the fluid are much smaller for turbulent flow than for 
laminar flow, as can be seen by comparing the "C M curves in Figures 57 and 66. 
The vertical distances between the minima and those points between the minima 
and maxima represent these differences. Due to the larger values of these 
differences for laminar flow, as well as the larger radial variations of fluid 
temperature, the temperature level in the mirror is generally considerably 
higher for laminar than for turbulent flow. 

The second figure in each triplet presents representative axial steady 
temperature distributions. The curve denoted X pertains to the column of 
nodes inside the shaded spot from node 3 at the top through node 155 at the 
bottom. The curve denoted Y pertains to the column of nodes just outside the 
shaded spot from node 5 at the top through node 153 at the bottom. The curve 
denoted Z pertains to the column of nodes centered at a radius of 8. 18 in. , 
from node 20 at the top to node 169 at the bottom. At this location, the thick- 
ness of the mirror is about 1. 4 in. 

The third figure in each triplet gives the variation with time of the nodes 
at the top and bottom of each of the columns for which results are presented in 
the second figure. The steady state temperatures of each of these nodes are 
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denoted by the locations of the corresponding symbols on the vertical axis. In 
general, the transient results are carried out to 0. 4 hours, the exception being 
those shown in Figure 59. 

Figures 57 through 59 present the results for the case of turbulent flow 
and a uniform fluid temperature. A uniform fluid temperature corresponds to 
an infinite flow rate, so this case is included primarily for its simplicity and 
for comparison purposes. It is observed from Figure 57 that the temperature 
difference between the front and back rows of nodes is roughly constant for 
radii greater than 3 inches. The drooping characteristic of curves F and B for 
radii approaching the periphery is due to shielding of the bevelled edge of the 
mirror. 

The constancy of the steady-state temperature difference across the 
present mirror is in contrast with the variable difference reported earlier for 
the tapered Cervit mirror, which was cooled at the back surface. (See, for 
example, Figure 32). This feature of the present results is apparently due to 
the uniform spacing of the cooling coils relative to the front surface of the 
mirror. It should also be noted that the average temperature difference 
between the front and back surfaces is considerably smaller with the present 
design than with the earlier one. The reason for the smaller temperature dif- 
ference becomes apparent upon consideration of the results in Figure 58. 

Curves Y and Z indicate that by far the major portion of the temperature differ- 
ence across the mirror occurs between the front surface and the coolant tubes, 
which are centered 1/2 inch below the front surface. This feature of the 
results is consistent with the geometry and the thermal boundary conditions. 

The negative temperature difference indicated by curve Z between the bottom of 
the tube at 0. 65 in. and the node in the bottom row at I. 2 in. is considerably 
larger than that obtained from the detailed analysis. (See, for example, Fig- 
ure 44,.) As a result, the temperature difference between the front and back 
rows of nodes is smaller in the present analysis than in the earlier one, This 
difference between the two sets of results is due to the coarser subdivisions 
used in the present analysis. The results of the detailed analysis are probably 
more accurate, so the present results may give optimistic estimates of the 
magnitudes of the thermal distortions. 

By comparing the results in Figures 57, 60, and 63, it may be noted that 
the temperature difference between the front and back rows of nodes at a 
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particular radiaL location is essentially independent of the radial variation in 
fluid temperature. This holds for laminar as well as turbulent flow, and it indi- 
cates that radial conduction is relatively unimportant. A similar feature was 
noted earlier in the results for the Cervit mirror. 

The transient results in Figures 59, 62, and 65 indicate that the tempera- 
tures at a radius of about 8. 2 in. are very near equilibrium 0. 15 hours after the 
start of illumination. Examination of the detailed results reveals that the ele- 
ments at smaller radii and near the front surface approach equilibrium almost 
as rapidly. However, the elements in the back row approach equilibrium much 
less rapidly. A feeling can be gained for the last statement by examining the 
curves labeled E in Figure 59, 62, and 65. In addition, the curves labeled A 
and D indicate that elements inside the shaded spot approach equilibrium slowly; 
and this is particularly the case for those elements at the bottom of the hub. 
Comparison of the results in Figures 59, 62, and 65 reveals that equilibrium is 
is approached at a rate, at a given node, which is essentially independent of the 
radial distribution of fluid temperature. 

Turning to the results for laminar flow in Figures 66, 67, and 68; 69, 70, 
and 71; and 72, 73, and 74, it is observed that the steady temperature difference 
across most of the mirror is about 0. 8°F larger than was the case for turbulent 
flow. This is due to the fact that the convective resistance is larger for laminar 
flow than for turbulent flow. As stated earlier, this is graphically illustrated 
by curve C in Figure 66. The curves in Figures 68, 71, and 74 show that equi- 
librium is approached much less rapidly with laminar flow than with turbulent 
flow. The reason is the same as above; the convective resistance is higher for 
the former than for the latter. 

A later design change came about in an effort to optimize the configuration 
with respect to distortion in the 1 -g field. This amounted to additional material 
1 in, thick at the hub and tapering to meet the original slope at a radius of 
9. 8 inches (Configuration No. 2, Figure 2). This new mirror configuration 
with the discretized elements and assigned node numbers 'is shown in Figure 75. 
These node numbers which are circled show the addition. The results of adding 
this material had a very slight effect on the temperature distribution. This 
can be seen by comparing Figures 80 and 84 which show the temperature pro- 
files of the front and back surfaces for configuration No. 1 and configuration 
No. 2 respectively. 
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Figures 76 through 87 show the front and back surface temperature 
profiles for those conditions which were analyzed for thermaL distortion. As 
indicated in Figure 76 the top curve is for the front surface and the bottom 
curve is for the back surface. 

Calculations have been performed to estimate the pressure drop of meth- 
anol water (80/20) through a circular cooling channel of diameter D and length L 
with constant radial spacing, Ar, between the centers of adjacent turns of a 
spiral pattern. The flow rates used are those required to maintain a 1°F or 2°F 
temperature difference between inlet and outlet of a 65-cm-diam. mirror with a 
12-cm-diam, shaded spot. The mirror absorbs 0. 14 of one solar constant. A 
standard equation for pressure drop was used, namely 


Ap 




(2) 


where f is the so-called Darcy-Weisbach friction factor. This equation does 
not account for curvature in the tube, but this is not thought to be important. 
The friction factor for a smooth tube was used. No allowance was made for 
pressure drops through fittings. 

The following equation was derived for the total length of the cooling 
channel for a mirror with inner radius r. and outer radius r^ and for a spacing 
between adjacent turns of r: 
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2-rr 
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The results are summarized below. 
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Tube 
Diameter 
(in. ) 

Radial 
Spacing 
(in. ) 

Total 

Channel 

Length 

(ft.) 

1/8 

3/8 

110 

3/16 

3/8 

no 

3/16 

9/16 

73. 4 

1/4 

3/8 

no 

1/4 

3/4 

55 

1/4 

3/4 

55 


Temp. 

Diff. 

( ° F) 

Mass Flow 
Rate 
(lb /in. 2 ) 

Pressure 
Drop 
(lb/in. 2 ) 

2 

143 

171 

2 

143 

24. 4 

2 

143 

16. 3 

2 

143 

6. 32 

2 

143 

3. 16 

1 

286 

10. 3 


Since turbulent flow would be a greater problem with respect to vibrations 
due to the fluid flow, it has been decided that laminar flow will be used. An 
estimated internal pressure of 16 psia was considered to be the maximum 
pressure internal to the mirror. 

This pressure has been based on the design where the pump is placed 
close to the outlet where the fluid exits the mirror. The minimum pressure of 
the system would be at this point, i. e. , the outlet of the mirror and the pump 
itself; therefore, the system must be designed so that the pressure at this point 
is above the boiling point of the coolant. 

A study was made to determine if it was feasible to use radiative cooling 
in order to effect the back surface temperatures which has been studied. Fig- 
ures 88 and 89 show the required sink temperature (T^) for the back surface of 
the mirror to be at 70, 80, 90 and 100°F, with a view factor between the mirror 
and sink varying from 0. 5 to 1. 0. The inset of Figure 88 shows a means of 
obtaining a view factor of 1. 0 between the mirror and the cold plate. The 
difference between the two models is that the solar absorptivity of Figure 88 is 
0. 14 and that for 89 is 0. 11. 

For the mirror to have a back surface temperature of 7 0°F (T^ of Fig- 
ure 88) with view factor of 1. 0, the sink plate temperature would have to be 
maintained at -40°F (T^,); therefore, a refrigeration system would be required. 

An investigation was made for a refrigeration system which could be 
used, but none was found available. The major problem with the system is 
that the compressor for the refrigeration cycle operates on the refrigerant 
vapor only. In a 1-g field the liquid can be separated from the vapor. In the 
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space application no practical method has been devised to ensure that only vapor 
goes into the compressor. 

For these reasons this method of controlling the mirror temperature is 
not feasible. 

Figure 16 shows, among other thingsi a temperature mapping of a steady - 
state computer analysis. This is the method utilized to transfer the tempera- 
ture data to the stress analysis group. These temperatures are from the 
printed listings on Table 14, 

Much of this written report has been taken directly from Dr, M.W. Wiidin's 
memos which he wrote while working on the primary mirror. 
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Figure 70. 
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Figure 71. 
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Figure 73. 
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Figure 76. Laminar Flow, Uniform 70° F Coolant Temperature, ULE Fused SiC> 2 Configuration Number 1 


TEMPERATURE, °F 




TEMPERATURE, °F 





/ 


Figure 79. Laminar Flow from Outside to Inside, ULE Fused Si0 2 , Configuration Number 1 
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Figure 80. Laminar Flow form Inside to Outside, ULE Fused SiO, 
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Figure 83. 


Temperature Gradient on Mirror Back 3°F Rise from Inside to Outside 
Beryllium Configuration Number 2 
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Figure 85. Laminar Flow from Inside to Outside, ULE Fused Si02> 
Configuration Number 2 
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Figure 87. Laminer Flow from Inside to Outside Mirror Front Face 0. 35 in. Thick at Mirror Center, 
Increasing Linerarly to 0. 50 in. Thick at the Mirror O. D. ULE Fused SiO^ Configuration 
Number 2 


Figure 88. 
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Figure 89. 
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